1/3 


AD-R189  323 
UNCLASSIFIED 


CORPUS  CHRISTI  INNER  HARBOR  SHOALING  INVESTIQATION(U> 
CAVER  (TROV  V)  DOVER  NJ  T  H  SMITH  ET  AL.  SEP  87 
MES/TR/HL-87-13 

F/G  13/2 


1.25  ■  1.4 


Destroy  this  report  when  no  longer  needed.  Do  no:  return 
it  to  the  originator. 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated 
by  other  authorized  documents. 


The  contents  of  this  report  are  not  to  be  used  for 
advertising,  publication,  c  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of 
such  commercial  products. 


Unclassified 


Hw.n.-nnt 


IFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


3  distribution  /availability  of  report 

Approved  for  public  release;  distribution 
unlimited . 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


la.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
Technical  Report  HL-87-13 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

USAEWES  (lf  'pp,,c*bW 

Hydraulics  Laboratory  WESHE-? 


6c  ADDRESS  [City.  State,  a nd  ZIP  Cod*) 

PO  Box  631 

Vicksburg,  MS  39180-0631 


7b  ADDRESS  (City.  State,  and  ZIP  Cod*) 


Ba  NAME  OE  FUNDING  /SPONSORING 
ORGANIZATION 

USAED,  Galveston 


8c  ADDRESS  (City,  Sttt*.  tnd  ZIP  Cod*) 

PO  Box  1229 

Galveston,  TX  77553-1229 


8b  OFFICE  SYM80L 
(If  *ppl«ablt) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 


PROJECT 

NO 


TASK 

WORK  UNIT 

NO 

ACCESSION  NO 

1 1  TITLE  (Intlud*  Security  Classification) 

Corpus  Christ!  Inner  Harbor  Shoaling  Invest igation 


12  PERSONAL  AUTHOR(S) 

Smith,  Tamsen  M. ;  McAnally,  William  H. ,  Jr.;  Teeter,  Allen  M. 


13a  TYPE  OF  REPORT  |13b  TIME  COVERED  |14  DATE  OF  REPORT  IYe*r,  Month,  D*y) 

Final  report  I  from _ to _  I  September  1987 


16  supplementary  notation 

Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road,  Springfield, 
VA  22161. 


17 _ CQSATI  COOES _ I  18  SUBJECT  TERMS  (Contmu*  on  reverie  if  n* ternary  and  identify  by  bJcxJr  number) 

Field  [  GROUP  |  sue  group  |  Corpus  Christ!  Bay  Numerical  modeling 

Corpus  Christ!  Harbor  Sedimentation 

Estuaries 


19  ABSTRACT  (Continue  on  reverie  if  neceuary  and  identify  by  block  number) 

A  combination  of  numerical  models  was  used  to  test  alternatives  for  shoaling  preven¬ 
tion  In  Corpus  Chrlstl  Harbor,  Texas.  The  vertically  averaged  model  system,  TABS-2,  was 
used  to  simulate  contributions  of  sediments  by  bay  waters  to  the  sediment  load.  The  lat¬ 
erally  averaged  estuarine  model,  I.AEMSED,  was  used  to  simulate  density  currents  In  the 
channel  and  sedimentation  that  occurs  at  the  harbor  entrance. 

Applications  of  the  models  testing  advance  maintenance,  removal  of  Industrial  dis¬ 
charges  and  withdrawals,  advance  maintenance  in  conjunction  with  a  sill,  and  movement  of 
the  disposal  areas  showed  a  20  percent  decrease  In  shoaling  as  a  result  of  Industrial 
activity  removal,  a  75  percent  decrease  In  sediments  entering  the  hay  channel  due  to  dis¬ 
posal  area  relocation,  and  practically  no  effect  on  shoaling  rates  resulting  from  advance 
maintenance . 

(font  f  nuedl 


2!  ABSTRACT  SECURITY  CLASSIFICATION 
Unc 1  ass  1  fled 


22b  telephone  (Induct*  Area  Code)  22c  OFFICE  SYMBOL 


20  distribution  >  availability  of  abstract 

Eunclassifieo/unumited  □  same  as  r pt 

□  DTiC  USERS 

22a  NAME  OF  RESPONSE  E  INDIVIDUAL 

I'nclasslf  led 

cwSEwSriS  «>  rw»  >«— 


19.  ABSTRACT  (Continued). 

Appendix  A  presents  the  results  of  a  reconnaissance  survey  on  shoaling  conditions  in 
Corpus  Chrlstl  Harbor.  Appendix  B  describes  the  TABS-2  numerical  modeling  system,  and 
Appendix  C  describes  the  theoretical  aspects  of  LAEMSED. 


_ Unc 1 assl f led _ 

iscuaiTv  classification  or  THIS  r»o* 


1*.  At.  iL  kL  itx  i^iL  i!  »1  *rt 


PREFACE 

This  report  describes  and  presents  the  results  from  the  Corpus  Christ! 
Harbor  Shoaling  Investigation  performed  by  the  US  Army  Engineer  Waterways  Ex¬ 
periment  Station  (WES)  for  the  US  Army  Engineer  District,  Galveston  (SWG) . 

The  study  was  conducted  by  personnel  of  the  Hydraulics  Laboratory  of  WES  under 
the  general  supervision  of  Mr.  F.  A.  Herrmann,  Jr.,  Chief  of  the  Hydraulics 
Laboratory,  Mr.  R.  A.  Sager,  Assistant  Chief  of  the  Hydraulics  Laboratory, 

Mr.  W.  H.  McAnally,  Chief  of  the  Estuaries  Division,  and  Mr.  G.  M.  Fisackerly, 
Chief  of  the  Estuarine  Processes  Branch  (EPB) .  Mr.  J.  M.  Keislich,  SWG,  was 
project  coordinator.  Work  on  the  project  began  in  December  1984  and  was 
completed  in  December  1986. 

Me  T.  M.  Smith,  EPB,  was  Project  Engineer.  Mr.  A.  M.  Teeter,  EPB, 
conducted  the  reconnaissance  survey  and  provided  consultation  for  the  field 
program  and  the  numerical  modeling.  Mr.  H.  A.  Benson,  EPB,  supervised  the 
field  data  collection  program,  and  Ms.  B.  P.  Donnell,  Estuarine  Simulation 
Branch  (ESB) ,  assisted  with  the  field  data  analysis.  Mr.  D.  P.  Bach,  ESB,  and 
Dr.  B.  H.  Johnson,  Math  Modeling  Group,  Hydraulic  Analysis  Division,  provided 
advice  and  assistance  with  the  numerical  modeling.  Ms.  Smith,  Mr.  McAnally, 
and  Mr.  Teeter  wrote  the  main  text  of  this  report.  Mr.  Teeter  wrote 
Appendix  A,  Dr.  Johnson  wrote  Appendix  C,  and  staff  of  the  Hydraulics  Labora¬ 
tory  wrote  Appendix  B.  Mrs.  Marsha  C.  Gay,  Information  Products  Division, 
Information  Technology  Laboratory,  edited  the  report. 

COL  Dwayne  G.  Lee,  CE,  Is  the  Commander  and  Director  of  WES. 

Dr.  Robert  W.  Whalin  is  the  Technical  Director. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement 
(metric)  units  as  follows: 

_ Multiply _ 

cubic  feet 
cubic  yards 
degrees  (angle) 
feet 

gallons  (US  liquid) 
inches 

knots  (international) 

miles  (US  nautical) 

miles  (US  statute) 

ounces  (US  fluid) 

pounds  (force) -second 
per  square  foot 

pounds  (mass) 


used  in  this  report  can 


_ _ 

0.02831685 

0.7645599 

0.01745329 

0.3048 

3.785412 

25.4 

0.514444 

1.852 

1.609344 

0.02957353 

47.88026 

0.4535924 


be  converted  to  SI 

To  Obtain 
cubic  metres 
cubic  metres 
radians 
metres 

cubic  decimetres 

millimetres 

metres  per  second 

kilometres 

kilometres 

cubic  decimetres 

pascals-second 

kilograms 


CORPUS  CHRIST I  INNER  HARBOR 
SHOALING  INVESTIGATION 


PART  I :  INTRODUCTION 

Background 

1.  Corpus  Christi  Bay  is  a  shallow,  wind-dominated  system  on  the  east 
coast  of  Texas  and  is  sheltered  from  the  open  waters  of  the  Gulf  of  Mexico  by 
a  low,  narrow  strip  of  land  known  as  Mustang  Island.  Corpus  Christi  Harbor, 
located  at  the  northwestern  corner  of  the  bay  (Figure  1),  experiences  high 
shoaling  rates  in  the  harbor  entrance  area.  A  planned  deepening  project  for 
the  harbor  is  expected  to  increase  shoaling  rates  and  volumes. 

2.  The  harbor  forms  a  part  of  the  Corpus  Christi  Ship  Channel,  a  deep- 
draft  project  providing  access  from  the  gulf  to  several  ports  within  the 
Corpus  Christi  Bay  system.  The  channel  extends  from  deep  water  in  the  gulf 
through  a  natural  inlet  at  Aransas  Pass,  then  across  2 i  miles*  of  open  water 
to  the  harbor.  The  navigable  channel  then  extends  an  additional  8.7  miles  via 
five  turning  basins  and  connecting  channels  to  a  terminating  point  at  Viola 
Turning  Basin  at  the  head  of  the  project.  The  project  has  been  completed  to  a 
depth  of  45  ft  mean  low  water,  with  the  exception  of  the  upper  portion  of  the 
harbor  which  has  an  average  depth  of  40  ft. 

3.  The  8.7-mile-long  portion  of  the  waterway  from  the  harbor  entrance 
to  the  Viola  Turning  Basin  (Figure  2)  is  bounded  by  land  on  three  sides.  The 
harbor  project,  isolating  the  waterway  from  Nueces  Bay  and  the  Nueces  River  to 
the  north,  was  constructed  by  disposing  of  dredged  material  to  create  a  pro¬ 
tective  embankment  along  the  north  side  of  the  channel.  Thus,  with  the  excep¬ 
tion  of  water  from  the  surface  drainage  and  the  effluents  from  municipal  and 
industrial  use,  all  water  In  the  inner  harbor  comes  from  the  adjacent  portions 
of  Corpus  Christi  Bay  and  the  navigation  channel  serving  the  harbor. 

4.  Historically,  the  first  2  miles  of  the  inner  harbor  channel  (from 
the  breakwaters  at  the  basin  entrance  through  the  basin  itself)  have  had 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  4. 


relatively  high  shoaling  rates  of  about  700,000  cu  yd  per  year.  The  deepening 
of  the  channel  to  its  45-ft  project  depth  beginning  in  1978  and  still  contin¬ 
uing  has  raised  concerns  that  the  shoaling  rate  in  this  area  will  increase. 
Concerns  have  also  been  expressed  regarding  the  fate  of  dredged  material  under 
disposal  plans  for  the  authorized  project.  Extensive  amounts  of  dredged  mate¬ 
rial  have  been  placed  in  bay  and  upland  disposal  areas.  According  to  the 
Corpus  Christi  Port  Authority  (Nueces  County  Navigation  Commission  1979)  ,  open 
water  disposal  of  sediments  dredged  from  the  harbor  is  no  longer  considered 
acceptable  by  regulatory  agencies,  and  land  disposal  sites  are  limited  in  size 
and  number.  Alternatives  to  reduce  the  shoaling  quantity  and  thereby  reduce 
the  disposal  requirements  are  needed.  Therefore,  the  processes  controlling 
sediment  deposition  in  the  inner  harbor  area  need  to  be  defined  so  that  the 
alternatives  can  be  evaluated. 

5.  The  Committee  on  Tidal  Hydraulics  (CTH) ,  Corps  of  Engineers,  US 

Army,  performed  a  study  of  this  sedimentation  problem  in  1964  (CTH  1965).  The 
primary  finding  of  that  study  was  the  identification  of  a  water  circulation 
pattern  similar  to  a  density  current  phenomenon  caused  by  one  or  both  of  the 
following:  (a)  the  large  industrial  discharges  to  and  withdrawals  from  the 

harbor  and  (b)  a  higher  mean  density  of  water  in  Corpus  Christi  Bay  than  in 
the  inner  harbor.  CTH  recommended  additional  comprehensive  field  studies  to 
determine  with  more  certainty  the  processes  controlling  the  water  circulation, 
sediment  transport,  and  deposition  patterns. 

Purpose 

6.  The  objectives  of  this  study  were  to  determine  the  sources  and 
causes  of  shoaling  occurring  in  the  inner  harbor  and  to  evaluate  alternative 
methods  for  reducing  the  high  volume  of  maintenance  dredging  required  near  the 
harbor  entrance.  In  addition,  the  schedule  for  the  overall  study  required  a 
rapid  completion  of  the  model  tests. 


Approach 

7.  This  investigation  was  conducted  in  three  phases.  Phase  I  was  a 
short-term  reconnaissance  survey  to  gather  the  hydraulic,  salinity,  tempera¬ 
ture,  and  sediment  data  necessary  to  characterize  the  system  and  to  identify 


phenomena  that  contribute  to  the  shoaling.  That  portion  of  the  study  was  com 
pleted  in  August  of  1984  and  appears  in  this  report  as  Appendix  A.  Phase  II 
activities  included  monitoring  harbor  and  bay  hydraulics,  salinities,  and  sus¬ 
pended  sediment  concentrations  under  a  variety  of  conditions  to  define  the 
factors  that  contribute  significantly  to  harbor  shoaling  and  to  support  the 
verification  of  numerical  models.  The  purpose  of  the  Phase  III  portion  of  the 
study  was  to  numerically  simulate  the  sedimentation  process  that  occurs  in  the 
harbor  to  evaluate  plans  for  the  reduction  of  harbor  shoaling. 

8.  Because  of  the  time  restriction,  limited  verification  of  the  numeri¬ 
cal  models  was  determined  to  be  sufficient  for  the  purposes  of  this  study. 

9.  A  combination  of  models  was  used  in  the  numerical  simulations.  The 
Phase  I  reconnaissance  survey  showed  that  the  predominant  shoaling  mechanism 
is  a  density-driven  current  along  the  bay  channel  that  brings  fine-grained 
sediments  into  the  harbor  entrance  area  where  the  sediments  deposit  and  ac¬ 
cumulate.  The  vertical  variation  in  velocities  required  a  vertical  model  (one 
that  reproduces  the  vertical  velocity  profiles  accurately)  to  simulate  sedi¬ 
mentation  processes  that  occurred.  LAEMSED  (Laterally  Averaged  Estuarine 
Model  with  Sediment)  was  used  to  model  the  ship  channel  and  harbor.  Modeling 
the  sediment  contributions  by  sources  in  the  bay  required  the  use  of  a  hori¬ 
zontal  model  (which  is  two-dimensional  in  the  horizontal  plane) .  The  numeri¬ 
cal  modeling  system,  TABS-2,  was  used.  Results  from  the  TABS-2  (bay)  model 
were  used  as  input  to  the  LAEMSED  (harbor)  model. 


PART  II:  DESCRIPTION  OF  MODELS 


Numerical  Modeling 

10.  Numerical  modeling  employs  special  computational  methods,  such  as 
iteration  and  approximation,  to  solve  mathematical  expressions  that  describe 
physical  phenomena.  Numerical  models  used  in  coastal  hydraulic  problems  are 
of  two  principal  types:  finite  difference  and  finite  element.  The  finite 
difference  method  approximates  derivatives  by  differences  in  the  values  of 
variables  over  finite  intervals  of  space  and  time.  This  method  requires  dis¬ 
cretization  of  space  and  time  into  regular  grids  of  computation  points.  The 
finite  element  method  employs  piecewise  approximations  of  mathematical  expres¬ 
sions  over  a  number  of  discrete  elements.  The  assemblage  of  piecewise  approx¬ 
imations  is  solved  as  a  set  of  simultaneous  equations  to  provide  answers  at 
points  in  space  (nodes)  and  time. 

11.  Numerical  models  are  classified  by  the  number  of  spatial  dimensions 
over  which  variables  are  permitted  to  change.  Thus  in  a  one-dimensional  flow 
model,  currents  are  averaged  over  two  dimensions  (usually  width  and  depth)  and 
vary  only  in  one  direction  (usually  longitudinally) .  Two-dimensional  models 
average  variables  over  one  spatial  dimension,  either  over  depth  (a  horizontal 
model)  or  width  (a  vertical  model).  Three-dimensional  models  solve  equations 
accounting  for  variation  of  the  variables  in  all  three  spatial  dimensions. 

12.  Numerical  models  are  capable  of  simulating  processes  that  cannot  be 
handled  in  any  other  way.  Once  a  numerical  model  has  been  formulated  and 
verified  for  a  given  area,  it  can  quickly  provide  results  for  different 
conditions . 


The  TABS-2  Modeling  System 

13.  The  TABS-2  system,  developed  by  the  Hydraulics  Laboratory,  US  Army 
Engineer  Waterways  Experiment  Station,  is  a  set  of  generalized  computer  pro¬ 
grams  for  two-dimensional  numerical  modeling  of  open-channel  flow,  transport 
processes,  and  sedimentation  problems  in  rivers,  reservoirs,  bays,  and  estu¬ 
aries  (Thomas  and  McAnally  1985a)  .  The  two  major  components  of  the  system 
used  in  this  study  were  RMA-2V,  the  finite  element  model  for  hydrodynamics, 
and  STUDH,  the  finite  element  sediment  model. 


14.  RMA-2V  solves  the  depth-integrated  equations  of  conservation  of 
mass  and  momentum  in  two  horizontal  directions.  Friction  is  calculated  with 
Manning's  equation  and  eddy  viscosity  coefficients  are  used  to  define  lateral 
turbulence  exchanges.  The  sedimentation  model,  STUDH,  solves  the  convection- 
diffusion  equation  with  bed  source  terms.  These  terms  are  structured  for 
either  sand  or  cohesive  sediments.  Clay  erosion  is  based  on  work  by  Par- 
theniades  (1962)  and  the  deposition  of  clay  utilizes  Krone's  (1962)  equations. 
The  hydrodynamics  from  RMA-2V  are  used  as  input  to  STUDH.  A  more  detailed 
description  of  both  models  appears  in  Appendix  B. 

The  LAEMSED  Model 

15.  LAEMSED  is  the  sediment  version  of  the  original  estuarine  model, 
LAEM  (Laterally  Averaged  Estuarine  Model) ,  developed  by  Edinger  and  Buchak 
(1981).  LAEM  was  developed  for  computing  stratified  flows  in  estuaries. 
Additions  to  LAEM  for  suspended  sediment  computations  and  to  include  a  sedi¬ 
ment  bed  model  were  made  by  Johnson  (Johnson,  Trawle,  and  Kee  in  preparation). 
These  modifications  were  patterned  after  work  by  Thomas  and  McAnally  (1985a) 
in  the  development  of  the  two-dimensional  sediment  model,  STUDH. 

16.  The  basic  flow  and  transport  equations  that  are  solved  in  LAEMSED 
are  statements  of  the  conservation  of  mass  and  momentum  of  the  flow  field  plus 
the  conservation  of  the  heat,  salt,  and  suspended  sediment  in  the  water  body. 
Boussinesq's  eddy  coefficient  concept  is  employed  to  account  for  effect  of 
turbulence  in  the  flow  field.  The  routines  that  compute  the  exchange  of  sedi¬ 
ment  between  the  sediment  bed  and  the  water  column  are  modifications  of  those 
found  in  STUDH.  A  more  detailed  description  of  the  model  is  given  in 
Appendix  C. 
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PART  III:  FIELD  INVESTIGATIONS 


17.  Field  data  obtained  in  Phases  I  and  II  of  the  study  were  used  for 
limited  model  verification.  The  data  collected  in  the  Phase  II  portion  con¬ 
sisted  of  long-term  tidal,  current,  and  water  sample  data  and  two  intensive 
surveys.  Long-term  data  collection  began  in  December  1984  and  was  completed 
in  October  1985.  Ten  months  of  tidal  data,  9  months  of  velocity  data, 

7  months  of  wind  data,  and  4  months  of  suspended  sediment  and  salinity  data 
were  collected.  Long-term  stations  are  shown  in  Figures  3  and  4.  In  addi¬ 
tion,  two  intensive  surveys  were  conducted  in  the  inner  harbor.  The  first, 
conducted  in  February  1985,  measured  velocities,  salinities,  and  suspended 
sediment  concentrations  at  nine  stations  over  a  25-hr  period.  The  second 
intensive  survey  was  conducted  in  June  1985.  Only  14  hr  of  data  were  col¬ 
lected  during  the  second  survey  because  of  hazardous  weather  conditions.  The 
field  data  collection  and  analysis  are  described  by  Smith  (in  preparation). 

18.  Because  the  current  velocities  measured  during  the  February  inten¬ 
sive  survey  were  below  the  threshold  of  the  current  meters  (approximately 
0.1  ft/sec)  a  majority  of  the  time,  only  the  data  from  the  June  survey  were 
used  for  model  verification.  Station  locations  for  this  survey  are  given  in 
Figure  5. 

19.  A  statistical  analysis  was  performed  on  wind  data  collected  at  the 
Corpus  Christi  Area  Office,  US  Army  Engineer  District,  Galveston,  and  on  cur¬ 
rent  velocity  data  measured  at  sta  V4  during  May  1985,  the  longest  period  for 
which  both  wind  and  velocity  data  near  the  harbor  entrance  were  available. 

Data  were  analyzed  to  determine  the  correlation  between  wind  velocities  and 
current  velocities  near  the  harbor  entrance.  Polar  histograms  of  velocities 
and  wind  speeds  are  shown  in  Figure  6.  Forty-one  percent  of  the  current 
velocities  occurring  at  this  station  had  direction  headings  in  the  north¬ 
western  quadrant  (Q4)  and  thirty-six  percent  of  the  headings  fell  in  the 
southwestern  quadrant  (Q3) .  Fifty-four  percent  of  the  wind  headings  were  in 
the  northwestern  quadrant  (from  the  southeast)  and  thirty-one  percent  occurred 
in  the  southwestern  quadrant.  Observations  of  the  data  from  other  times  of 
the  year  show  the  same  pattern.  These  data  and  the  physical  characteristics 
of  the  system  (small  tide  range  and  wide  shallow  bay)  suggest  that  currents 
near  the  channel  at  sta  V4  are  strongly  influenced  by  wind. 

20.  It  is  reported  that  during  most  of  the  year,  winds  at  Corpus 
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Christi  are  generally  from  the  southeast  (Texas  A&M  Research  Foundation  1973) . 
Also,  throughout  the  period  of  record  (January-October  1985),  a  majority  of 
the  current  velocities  were  toward  the  northwest.  Therefore,  it  is  reasonable 
to  assume  that  a  majority  of  the  time,  current  directions  in  the  bay  near  the 
channel  are  heavily  influenced  by  winds  and  flow  is  predominantly  toward  the 
entrance . 

21.  Present  dredging  practices  at  the  Corpus  Christi  Ship  Channel  were 
defined  during  field  data  collection.  Material  that  is  dredged  in  the  bay 
channel  is  discharged  in  open  water  about  1,600  ft  south  of  the  channel.  The 
mound  from  these  point  discharges  has  never  been  observed  to  reach  the  sur¬ 
face.  This  observation  confirms  a  prediction  stated  in  the  Phase  1  report 
(Appendix  A,  paragraph  47)  that  the  dredged  material  is  spreading  over  a  large 
surface  area  before  mounding,  becoming  easily  erodible. 

22.  Observations  of  suspended  sediment  concentrations  measured  from 
automatically  withdrawn  water  samples  (locations  shown  in  Figure  4)  show  that 
in  the  earlier  months  of  data  collection  (June-July  1985)  the  highest  concen¬ 
trations  measured  occurred  near  sta  SS4  on  the  southern  side  of  the  channel. 
Concentrations  at  this  station  declined  during  the  latter  portion  of  the  sur¬ 
vey  period.  The  probable  reason  for  the  elevated  values  in  suspended  sediment 
concentrations  is  the  dredging  of  the  bay  channel  taking  place  near  the  sta¬ 
tion  location  at  that  time.  Concentrations  at  sta  SSI  near  Corpus  Christi 
Beach  remained  relatively  high  (about  120  ppm)  during  the  survey  period. 

23.  Hydrographic  surveys  from  1965  to  1984  were  obtained  from  the 
Corpus  Christi  Area  Office.  Shoaling  rates  were  computed  by  comparing  the 
postdredge  survey  with  the  following  predredge  survey,  thus  establishing  the 
amount  of  material  that  had  shoaled  during  that  period.  The  following  tabula¬ 
tion  shows  the  results  of  that  comparison.  Shoaling  rates  from  0.7  to  2.0  ft 
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per  year  were  computed  for  the  harbor  entrance  area  and  the  Corpus  Christ i 
Turning  Basin.  The  hydrographic  surveys  near  the  disposal  areas  did  not  re¬ 
veal  any  significant  accumulation  of  mounded  sediment. 

24.  Corpus  Christi  Bay  supports  a  productive  shrimp  fishery,  and  the 
use  of  weighted  nets  may  make  a  significant  contribution  to  the  amount  of  sus¬ 
pended  solids.  Schubel  et  al.  (1978)  report  that  the  amount  of  sediment  re¬ 
suspended  each  year  by  shrimpers  is  approximately  16  to  131  times  that  dredged 
for  channel  maintenance.  During  the  reconnaissance  survey  and  field  data 
collection  efforts  it  was  noted  that  shrimping  was  particularly  heavy  near  the 
harbor  entrance. 
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PART  IV:  VERIFICATION  OF  THE  MODELS 


RMA-2V  Verification 


25.  As  previously  stated,  verification  of  the  models  was  limited  due  to 
time  constraints.  The  finite  element  mesh  of  the  bay  used  as  input  to  RMA-2V 
and  STUDH  is  shown  in  Figure  7.  The  mesh  was  made  up  of  elements  (three-  or 
four-sided  figures)  and  nodes  (corners  and  midside  points  of  the  elements) . 

The  final  numerical  mesh,  consisting  of  892  elements  and  2,810  nodes,  con¬ 
tained  a  detailed  representation  of  the  harbor  entrance  and  the  bay  disposal 
areas.  Elements  of  the  various  regions  of  the  bay  were  assigned  different 
types  for  specifying  hydrodynamic  and  sediment  parameters  (such  as  Manning's 

n  and  diffusion  coefficients).  This  network  contained  34  different  element 
types.  Boundary  conditions  for  the  hydrodynamic  model  consisted  of  velocity 
specifications  at  the  Nueces  River  (the  primary  source  of  fresh  water  to  the 
system)  and  water-surface  elevations  at  the  ocean  boundary. 

26.  RMA-2V  verification  was  performed  using  data  collected  at  the  sta¬ 
tions  in  the  bay  during  15-16  June  1985.  Tidal  data  were  filtered  to  remove 
meteorological  effects  and  high-frequency  oscillations.  The  mean  water- 
surface  elevation  was  obtained  by  filtering  the  9  months  of  available  tidal 
records.  This  elevation  was  assigned  a  value  of  100.0  in  RMA-2V  so  that  all 
bed  elevations  and  water-surface  elevations  would  be  positive.  The  tide  range 
at  the  exit  (Port  Aransas)  during  this  period  was  1.2  ft.  The  Nueces  River 
average  discharge  during  June  was  620  cfs.  Hydrodynamic  verification  of 
RMA-2V  was  accomplished  by  adjusting  bottom  roughness  and  eddy  viscosity  coef¬ 
ficients,  then  comparing  model  water-surface  elevations  and  velocities  with 


prototype  data.  One-hour  time-steps  were  used.  Eddy  viscosity  coefficients 

2  2 


ranged  from  10  lb-sec/ft  for  the  very  small  elements  to  400  lb-sec/ft^  for 
the  extremely  large  elements.  The  following  tabulation  provides  a  summary  of 
the  selected  verification  roughness  values  assigned  to  various  regions: 


Region 


Manning' s  n 


Ship  Channel 
Corpus  Christi  Bay 
Nueces  Bav 


0.012 

0.025-0.04 

0.06 


After  sufficient  time  was  allowed  for  mode  1  spin-up,  a  26-hr  test 


was  conducted.  Plates  1-3  show  prototype-numerical  model  wat or -sur 1  are 
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elevations  at  the  six  tide  stations,  and  Plates  4-6  show  prototype-numerical 
model  velocity  at  three  velocity  stations.  Phase  and  range  agreement  were 
good  at  the  tide  stations  near  the  area  of  interest  (sta  TG1-TG4) .  With  the 
exception  of  a  few  data  points,  good  agreement  was  also  obtained  at  the  ve¬ 
locity  stations.  At  lower  velocities  (those  below  0.1  ft/sec)  the  magnitude 
of  the  current  was  below  the  threshold  of  the  meter,  and  the  field  data  tended 
to  be  inconsistent  at  some  points. 

28.  After  adjustment  of  the  bay  hydrodynamic  model  to  the  15-16  June 
data,  a  test  was  conducted  corresponding  to  the  time  of  the  harbor  survey 
(17-18  June).  The  tide  range  at  Port  Aransas  during  this  period  was  1.3  ft. 
Plates  7-9  show  prototype-numerical  model  water-surface  elevation  comparisons 
and  Plate  10  shows  prototype-numerical  model  velocity  comparisons  for  stations 
where  field  data  were  available.  Plots  showing  hourly  current  velocities 
throughout  the  tidal  cycle  are  shown  in  Plates  11-17. 

29.  As  stated  in  paragraph  19,  bay  currents  are  strongly  dependent  on 
wind,  especially  during  periods  of  low  tidal  energy.  The  verification  period 
was  during  a  period  of  strong  tidal  energy,  and  running  RMA-2V  with  wind  did 
not  improve  the  verification  over  the  no-wind  case.  Resources  and  time  for 
the  feasibility  study  did  not  allow  running  the  various  seasonal  wind  condi¬ 
tions  prevalent  in  the  area;  therefore,  wind  was  not  considered  in  the  plan 
runs.  Model  runs  with  varying  wind  speeds  and  directions  should  be  included 
in  the  design  study. 

STUDH  Verification 

30.  Shoaling  problems  in  the  harbor  result  primarily  from  the  deposi¬ 
tion  of  fine-grained  material.  Thus  the  sediment  was  considered  to  be  a  clay. 
Because  of  the  required  study  schedule,  a  very  limited  verification  of  STUDH 
was  conducted.  Since  the  objective  of  this  portion  of  the  study  was  to  exam¬ 
ine  dredged  material  return  to  the  channel,  the  sediment  in  the  present  dis¬ 
posal  areas  was  allowed  to  erode  while  the  remainder  of  the  bay  was  assigned 
an  inerodible  bottom.  Parameters  such  as  critical  shear  stress  for  erosion, 
critical  shear  stress  for  deposition,  and  diffusion  coefficients  were  varied 
until  model  concentrations  at  the  three  bay  sediment  stations  matched  those 
observed  in  the  prototype.  Table  1  provides  the  values  used  for  the  various 
sediment  coefficients  in  the  model. 


31.  During  the  intensive  survey,  suspended  sediment  concentrations  re¬ 
mained  relatively  constant  over  the  tidal  cycle  with  the  exception  of  the  last 
few  hours  (storm  event) .  The  concentrations  computed  in  STUDH  at  nodes  cor¬ 
responding  to  these  prototype  stations  also  remained  fairly  constant  during 
the  last  half  of  the  tidal  cycle.  These  average  concentrations  for  model  and 
prototype  at  the  four  sediment  stations  are  given  in  the  following  tabulation. 

Prototype  STUDH 

Concentrations  Concentrations 


Station 

ppm 

ppm 

SSI 

51 

91 

SS2 

33 

35 

SS3 

16 

30 

SS4 

83 

73 

At  sta  SSI  near  Corpus  Christ!  Beach,  STUDH  concentrations  exceeded  those  mea¬ 
sured  in  the  field,  but  are  still  considered  reasonable.  Agreement  was  good 
at  the  other  three  stations. 

LAEMSED  Verification 

Hydrodynamic  verification 

32.  LAEMSED  modeled  the  Corpus  Christi  Ship  Channel  from  Viola  Basin  to 
approximately  8  miles  east  of  the  harbor  entrance  (near  Port  Ingleside) .  The 
grid  boundaries  from  the  Viola  Basin  to  the  harbor  entrance  were  the  land 
boundaries;  and  from  the  harbor  entrance  east  to  its  right  boundary,  the  grid 
encompassed  only  the  ship  channel.  A  schematic  of  a  portion  of  the  finite 
difference  grid  is  shown  in  Figure  8.  Associated  with  each  grid  cell  is  a 
width  corresponding  to  the  width  of  the  channel  in  that  area.  Each  cell  had  a 
length  dx  of  1,000  ft  and  a  height  dh  of  3  ft.  The  entire  grid  contained 
102  reaches  (dx's)  and  21  vertical  layers  (dh's).  One-minute  time-steps  were 
required . 

33.  Lateral  inflows  into  the  channel  from  the  bay  were  calculated  from 
RMA-2V  results.  The  flow  across  channel  element  sides  was  calculated,  and 
these  values  served  as  input  into  LAEMSED.  Associated  with  each  inflow  is  a 
sediment  concentration  obtained  from  field  data.  Industrial  inflows  QIN  and 
withdrawals  QWD  were  also  modeled  in  LAEMSED.  Information  concerning  two 
withdrawals  and  four  discharges  from  industries  along  the  harbor  were  used  as 

22 


■  f -r »  *'•  <,  «  -  "r<1  ^  .  .• ,  /  •  ,  ,* 


CORPUS  CHRIST! 


input  into  the  model.  These  are  given  in  Table  2.  Their  locations  are  shown 
in  Figure  2. 

34.  The  flow  computations  in  LAEMSED  were  verified  through  use  of  the 
field  data  collected  18  June  1985  and  discussed  in  Part  III.  The  tide  range 
at  the  right  boundary  (gage  TG5)  for  this  period  was  0.9  ft.  Salinities 
varied  during  the  observation  from  33.5  to  36.8  ppt  near  the  bottom  and  28.2 
to  33.5  ppt  at  the  surface.  A  generally  linear  vertical  distribution  of 
salinity  had  been  observed  in  the  field  and  was  used  in  the  model.  Sediment 
concentrations  varied  vertically  from  10  g /l  near  the  surface  to  30  g/ 2,  near 
the  bottom  but  did  not  vary  significantly  in  time  during  the  data  collection. 
Recorded  tides  at  five  locations  and  velocity  measurements  at  surface,  mid¬ 
depth,  and  bottom  at  six  locations  were  used  for  adjustment  purposes.  The 
eastern  boundary  of  the  grid  extended  across  the  bay  to  sta  TG5  where  water 
levels  were  recorded,  providing  water  level  boundary  conditions  for  the  model. 
Four  tidal  cycles  were  run  as  a  spin-up  period  before  computed  results  were 
compared  with  the  field  data.  A  comparison  of  computed  and  recorded  tides  at 
the  other  four  tide  stations  is  presented  in  Plates  18  and  19.  Phase  and 
range  comparison  at  all  stations  was  good  for  the  first  15  hr.  The  last  10  hr 
corresponds  to  the  storm  that  occurred  on  18  June.  The  prevalent  south¬ 
easterly  wind  shifted  to  a  northern  wind,  causing  the  depression  in  water  sur¬ 
face  observed  in  the  prototype  data.  The  storm  was  not  modeled. 

35.  Plots  of  computed  and  recorded  velocities  at  surface,  middepth,  and 

bottom  are  shown  in  Plates  20-25.  At  most  stations,  both  model  and  prototype 
data  show  a  net  flow  into  the  harbor  at  the  bottom  (flood  currents)  and  ouc  of 
the  harbor  at  the  surface  (ebb  currents) .  Stations  near  the  entrance  show  a 

net  inflow  at  middepth  as  well.  Phase  agreement  at  most  stations  was  gener¬ 

ally  good;  however,  prototype  magnitudes  were  much  greater  than  model,  espe¬ 
cially  at  stations  near  the  harbor  entrance.  As  noted  in  Appendix  A,  para¬ 
graph  39,  these  field  data  were  so  erratic  as  to  be  unreasonable.  High  winds 
during  the  survey  period  caused  motion  of  the  boats  from  whic’  measurements 

were  taken  and  may  have  induced  short-period  surges  that  contaminated  the 

data.  Velocities  measured  during  the  intensive  survey  in  February  were  more 
on  the  order  of  those  computed  by  the  model.  For  these  reasons  and  others, 
the  discrepancies  are  not  considered  to  reflect  a  poor  verification. 

Sediment  verification 

36.  As  in  STUDH,  only  the  sediment  parameters  for  clays  were  used  in 


LAEMSED.  Table  1  provides  the  values  used  for  the  various  sediment  coeffi¬ 
cients  in  the  model.  The  bed  was  assumed  to  be  initially  free  of  sediment. 
Initial  sediment  concentrations  in  each  layer  were  interpolated  from  field 
data  (paragraph  34) .  Sediment  was  introduced  into  the  model  by  assigning  a 
concentration  to  the  lateral  inflows  that  represent  flow  into  the  channel  from 
the  bay.  These  concentrations  were  also  obtained  from  field  data  (para¬ 
graph  33) . 

37.  The  shoaling  rates  computed  from  hydrographic  surveys  (para¬ 
graph  21)  were  for  two  separate  areas.  Region  1  is  the  4 ,000-f t-long  section 
of  channel  just  east  of  the  entrance,  and  region  2  consists  of  the  harbor 
entrance  and  the  Corpus  Christi  Turning  Basin.  The  tabulation  in  paragraph  23 
shows  the  computed  shoaling  rates  in  these  areas  compared  with  those  estimated 
from  the  hydrographic  surveys.  West  of  region  2,  further  into  the  harbor, 
both  prototype  and  model  shoaling  rates  were  less  than  0.25  ft  per  year.  The 
only  major  discrepancy  is  in  region  2  where  the  model  value  was  close  to  the 
rate  computed  during  the  first  period,  but  less  than  half  of  that  of  the  most 
recent  shoaling  rates.  During  this  period  of  time  (1980  to  1983),  more  than 
the  usual  amount  of  dredging  was  being  conducted  in  the  bay  channel.  Many 
sections  of  the  channel  were  being  dredged  to  the  new  authorized  depth  of 
45  ft.  This  activity  generated  more  dredged  material,  which  was  placed  in  the 
open  water  disposal  areas  and  could  have  gradually  flowed  back  into  the  chan¬ 
nel  and  from  there  into  the  harbor. 
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PART  V:  APPLICATION  OF  THE  MODELS 

38.  A  number  of  alternatives  were  identified  that  could  reduce  the 
annual  maintenance  dredging.  Presented  in  the  following  paragraphs  is  a 
description  of  each  of  the  tested  alternatives  and  how  the  two  models  were 
used  to  evaluate  their  effectiveness  in  reducing  maintenance  dredging. 

Modification  of  Industrial  Discharges  and  Withdrawals 

39.  One  of  the  factors  believed  to  be  contributing  to  the  density  cur¬ 
rent  phenomenon  observed  in  the  inner  harbor  is  the  large  amount  of  industrial 
discharges  into  and  withdrawals  from  the  project.  Therefore,  the  effect  of 
their  removal  on  the  circulation  pattern  and  shoaling  rate  was  tested.  Infor¬ 
mation  concerning  two  withdrawals  and  four  discharges  from  industries  along 
the  harbor  is  given  in  Table  2.  Their  locations  are  shown  in  Figure.  2. 

AO.  The  Central  Power  and  Light  Company  (Nueces  Plant)  withdraws  about 
A50  mgd  for  cooling  purposes  from  the  inner  harbor  and  discharges  this  water 
into  Nueces  Bay.  This  diversion  of  water  represents  a  sink  within  the  harbor 
and  must  be  replaced  with  waters  from  Corpus  Christi  Bay.  Because  of  the  mix¬ 
ing  mechanisms  and  vertical  density  gradients  at  the  harbor  entrance,  it  is 
possible  that  the  waters  that  are  brought  into  the  harbor  to  replenish  the 
cooling  water  diversions  are  the  saltier,  denser,  sediment-laden  waters  near 
the  bottom  of  the  channel.  Modification  of  the  cooling  water  diversion, 
either  through  a  reduction  of  the  volume  removed  or  through  discharging  the 
cooling  water  back  into  the  harbor  following  its  use,  would  reduce  or  elimi¬ 
nate  the  need  for  bay  waters  to  make  up  for  the  diversion. 

41.  In  addition  to  the  large  withdrawal  from  the  channel,  historical 
records  show  that  an  average  daily  freshwater  discharge  of  about  18  mgd  is 
released  into  the  upper  harbor.  The  mechanisms  of  freshwater/saltwater  mixing 
could  be  such  that  the  freshwater  discharges  induce  a  flow  into  the  harbor 
entrance  from  the  saltier,  denser,  sediment-laden  waters  of  the  bay  channel. 

As  the  bay  channel  waters  enter  the  harbor  entrance,  the  sediment  particles 
could  settle  to  form  the  harbor  entrance  shoals.  These  inflows  and  with¬ 
drawals  were  simulated  in  the  LAEMSF.D  model,  and  their  effects  on  the  circula¬ 
tion  pattern  in  the  harbor  and  on  the  sediment  deposition  were  tested. 
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Advance  Maintenance  Dredgins 


42.  Advance  maintenance  dredging  is  the  practice  of  dredging  channel 
dimensions  beyond  project  requirements  to  provide  additional  area  within  the 
project  limits  for  sediments  to  deposit  and  increase  the  time  between  dredg¬ 
ing.  If  the  additional  channel  dimensions  result  in  reducing  the  costs  and 
frequency  of  maintenance  dredging,  then  this  option  is  economically  feasible. 
Using  LAEMSED,  the  change  in  shoaling  rates,  and  therefore  the  change  in  the 
required  dredging  frequency,  as  a  result  of  the  advance  maintenance  was 
tested.  This  alternative  would  not  achieve  the  objective  of  reducing  the 
total  volume  dredged,  but  could  reduce  dredging  costs  and  provide  flexibility 
in  disposal. 

Structural  Modifications 

43.  Structural  modifications  near  the  harbor  entrance  is  a  possible 
alternative  for  diverting  sediment-laden  waters  from  the  harbor  entrance.  A 
sill  near  the  channel  bottom  might  act  to  inhibit  encroachment  of  saline  and 
sediment-laden  waters  into  the  harbor  along  the  channel  bottom.  Simulations 
of  such  a  structure  near  the  channel  bottom  in  conjunction  with  advance  main¬ 
tenance  were  performed,  using  LAEMSED,  to  observe  the  effect  that  the  restric¬ 
tion  of  the  denser  bottom  layers  would  have  on  sedimentation. 

Relocation  of  Dredged  Material  Disposal  Areas 


44.  Disposal  practices  at  Corpus  Christi  Harbor,  discussed  In  Part  Til, 
are  believed  to  contribute  directly  to  the  shoaling  problem  that  occurs  there. 
As  mentioned  in  paragraph  23,  surveys  conducted  at  the  disposal  areas  show 
little  evidence  of  dredged  material  remaining  at  the  disposal  areas.  Evidence 
was  observed  in  the  reconnaissance  survey  that  these  materials  are  possiblv 
being  returned  to  the  channel  by  wind-generated  currents  and  bv  the  general 
circulation  of  the  bay  (Appendix  A,  paragraphs  50  and  53).  Currents  in  shal¬ 
low  areas  adjacent  to  the  channel  flow  obliquely  to  the  channel.  As  the 
sediments  are  introduced  into  the  channel  by  this  oblique  current,  densitv 
mechanisms  can  move  them  to  the  harbor  entrance  where  thev  redeposit  to  form 
the  recurring  shoaling.  Using  the  TABS-2  model,  these  sites  were  relocated  to 
the  north  and  changes  in  the  amount  of  sediment  entering  the  channel  were 


PART  VI:  PLAN  TESTING 


Conditions  Tested 


45.  To  evaluate  the  impact  of  the  alternatives  discussed  in  Part  V, 
five  plan  tests  were  performed.  They  are  as  follows: 

a.  Plan  1.  Elimination  of  industrial  withdrawals  and  discharges 
(LAEMSED) . 

b.  Plan  2.  Advance  maintenance  dredging  (LAEMSED). 

c.  Plan  3.  Advance  maintenance  dredging  in  combination  with  a 
sill  (LAEMSED) . 

d.  Plan  4.  Relocation  of  bay  disposal  sites  to  1,900  ft  north  of 
the  channel  (STUDH) . 

e.  Plan  5.  Relocation  of  bay  disposal  sites  to  6,000  ft  north  of 
the  channel  (STUDH) . 

Boundary  and  initial  conditions  for  plan  tests  were  identical  to  those  used  in 
the  verification. 

Plan  1  Results 

46.  Locations  of  the  cells  for  which  sedimentation  rates  are  shown  are 
given  in  Figure  9.  The  shoaling  rates  from  the  LAEMSED  base  condition  test  of 
17-18  June  and  those  from  a  simulation  with  alternative  1  modifications 
(described  in  paragraphs  39-41)  are  compared  in  Figure  10.  As  can  be  seen 


Figure  9.  Location  of  LAEMSED  computation  cells  near  harbor  entrance 


d.  Cell  50 

Figure  10.  (Concluded) 


from  these  graphs.  Flan  1  decreased  the  shoaling  rate  at  the  harbor  entrance 
and  In  the  bav  channel  just  east  of  the  entrance.  lust  west  of  t  fie  entrance 
the  shoaling  staved  about  the  same,  and  in  the  Corpus  Christ i  Turning  Basin, 
shoaling  rates  increased. 

47.  The  major  difference  Flan  1  made  in  the  velocities  was  in  t  fie  top 
half  of  the  channel  (since  most  of  the  discharges  and  withdrawals  were  at  the 
surface).  Base  and  Plan  1  surface  velocities  for  sta  A  and  !'  are  shown  in 
Figure  11.  The  flood  predominance  observed  near  the  bottom  at  stations  near 
the  entrance  in  the  base  conditions  was  not  affected.  In  the  eastern  end  of 
the  harbor  and  in  the  bav  channel,  the  surface  velocities  became  more  ebb- 
dominated  as  seen  in  Figure  11a.  Therefore,  less  sediment  was  drawn  into  the 
channel  and  shoaling  decreased.  However,  in  the  western  half  of  the  first 
turning  basin  and  further  into  the  harbor,  the  surface  velocities  shifted  more 
in  the  flood  direction  (see  sta  0 ,  Figure  lib).  The  removal  of  the  large  sur¬ 
face  withdrawal  resulted  in  a  reduction  of  vertical  velocities,  causing  parti¬ 
cle  settling  velocities  to  become  less  impeded.  Slightlv  more  sediment  was 
allowed  to  deposit,  leading  to  the  small  increase  in  shoaling  observed  in  that 
area.  From  the  western  end  of  the  turning  basin  to  the  breakwaters,  the  net 
change  was  a  decrease  in  total  shoaling  volume  of  approx imat e 1 v  percent. 

Plan  2  Results 

48.  Oiannel  bottom  elevations  in  the  l.AFMSFb  grid  were  deepened  bv  h  ft 
to  reflect  the  result  of  advance  maintenance  dredging.  Figures  17  and  1  1  show 
the  resulting  change  In  shoaling  rates  at  cells  located  at  the  area  of  inter¬ 
est.  Throughout  the  entrance  and  turning  basin,  shoaling  rate1-  increased 
slightlv.  The  maximum  increase  observed  was  almost  1  in.  per  vear,  and  t  h.  1 
was  in  the  turning  basin.  However,  the  bav  channel  east  of  the  entrain  e, 
where  base  shoaling  rates  were  the  highest,  experienced  the  smallest  amount  of 
increase  (about  0.1  in.  per  vear).  An  increase  of  less  than  1  in.  pet  vear  in 
an  area  where  the  base  condition  shoaling  rates  are  already  ^  in.  per  vear  m,o 
he  small  enough  to  make  this  plan  economically  attractive,  and  an  ei  n'ornfi 
analysis  is  appropriate.  If  t  fie  cost  reduction  can  be  obtained,  the  mo"ev 
saved  could  be  used  to  revise  disposal  operations,  as  discussed  late?. 
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Kivure  13.  l.AEMSED  seilimentat  Ion  results,  base  to  Plan  ?  anil  Plan  1 


Plan  3  Results 

49.  In  an  attempt  to  offset  the  increased  shoaling  experienced  in  the 
Plan  2  test,  a  3-ft-high  sill  was  added  in  conjunction  with  the  advance  main¬ 
tenance.  The  sill  was  added  at  the  channel  bottom  just  west  of  the  entrance 
(cell  51).  Since  the  sill  was  added  in  combination  with  advance  maintenance 
dredging,  a  50-ft  controlling  depth  at  mean  low  water  was  maintained.  The 
resulting  shoaling  rates  at  cells  49  and  50  are  shown  in  Figure  13.  Just  east 
of  the  sill,  shoaling  decreased  from  the  Plan  2  alternative  by  almost  1.5  in. 


per  year.  Beyond  this  point  further  into  the  harbor,  the  effect  became 
negligible . 


Plan  4  Results 


50.  Results  from  RMA-2V  showed  that  currents  generally  flowed  from  the 
existing  disposal  areas  into  the  channel.  Figure  14  shows  peak  ebb  and  flood 
velocities  in  the  northwestern  section  of  the  bay.  STUDH  was  used  to  observe 
the  movement  of  sediment  which  was  allowed  to  erode  from  the  existing  disposal 
areas  as  described  in  paragraph  30.  The  mass  of  eroded  material  flowing  into 
the  channel  area  from  the  disposal  area  was  calculated  by  multiplying  flow  by 
sediment  concentration  at  each  time-step.  The  average  flux  of  sediment  into 
the  channel  was  35,000  kg/hr.  Since  LAEMSED  results  showed  that  much  of  the 
material  entering  the  bay  channel  Is  pulled  into  the  inner  harbor,  a  means  of 
reducing  the  amount  of  material  entering  the  channel  is  needed.  Based  on  ob¬ 
servations  of  the  circulation  patterns  and  field  sediment  data,  3  probable 
solution  seemed  to  be  to  relocate  the  disposal  sites  to  the  north  of  the  chan¬ 
nel  at  the  same  distance  from  the  channel  (2,000  ft).  The  finite  element  mesh 
was  revised  to  reflect  the  new  disposal  areas  in  the  model .  The  same  sediment 
characteristics  used  in  the  base  test  were  assigned  to  the  plan  disposal  areas 
and  a  one-tidal-cyc.le  simulation  was  made.  Figure  15  shows  base  and  plan  dis¬ 
posal  site  locations.  A  graph  of  suspended  sediment  concentrations  in  the 
channel  at  sta  R1  (shown  in  Figure  15)  for  base  and  Plan  4  conditions  is  shown 
in  Figure  16.  The  concentrations  computed  in  the  Plan  4  simulation  reached 
almost  the  same  level  as  thev  did  in  the  base  during  the  ebb  phase  of  the 
tide.  Close  inspection  of  the  velocities  in  Figure  14  show  that  during  ebb, 
the  direction  of  flow  is  slightly  toward  the  channel  in  the  vicinitv  of  the 
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Figure  15.  Base  and  Plan  4  disposal  site  locations 
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BASE  CONDITION 


MOOEL  TIME,  HR 

Figure  16.  Sedimentation  results,  base  to  Plan  4,  sta  R1 

plan  disposal  areas.  Also,  velocities  are  slightly  larger  north  of  the  chan¬ 
nel  so  more  material  is  eroded.  Therefore,  there  was  not  a  significant 
decrease  (less  than  5  percent)  in  the  amount  of  sediment  flowing  into  the 
channel . 

Plan  5  Results 

51.  Due  to  the  results  of  the  Plan  4  tests,  it  was  decided  to  relocate 
the  plan  disposal  areas  further  north  for  additional  tests.  The  Plan  5  dis¬ 
posal  sites  are  shown  in  Figure  17.  Sites  further  north  were  chosen  because 
RMA-2V  results  showed  flow  in  this  area  to  be  almost  always  either  parallel  to 
or  away  from  the  channel.  Their  distance  from  the  channel  was  limited  to 
6,000  ft  for  dredging  cost  purposes.  Figure  18  shows  Plan  5  disposal  sites 
with  RMA-2V  velocity  fields.  As  can  be  seen,  the  current  velocities  at  the 
most  western  site  (site  1)  are  much  higher  than  at  the  other  three  sites  due 
to  the  influence  of  the  Nueces  Bay  tidal  prism.  A  6-hr  STUDH  simulation 
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showed  that  erosion  of  Plan  5  disposal  site  1  was  almost  twice  as  much  as  that 
of  any  of  the  other  base  or  plan  sites.  This  created  large  suspended  sediment 
concentrations  which,  even  before  the  ebb  phase  of  the  tide,  had  begun  to 
enter  the  channel  as  a  result  of  diffusion.  Therefore,  Plan  5  was  revised  so 
that  no  material  was  placed  in  site  1  and  more  material  was  placed  in 
sites  2-4.  This  was  then  labeled  as  Plan  5b.  A  25-hr  STUDH  test  allowing 
erosion  of  the  Plan  5b  disposal  areas  was  conducted.  Model  concentrations  at 
sta  R1  after  25  hr  were  negligible.  Figure  19  shows  STUDH  base  and  Plan  5b 
concentrations  along  the  channel  at  hour  25.  There  was  a  resulting  decrease 
of  75  percent  sediment  mass  flux  into  the  channel. 
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PART  VII:  CONCLUSIONS  AND  RECOMMENDATIONS 


52.  Analysis  of  pre-  and  post-dredging  surveys  indicate  that  annual 
shoaling  rates  in  the  inner  harbor  area  have  been  between  0. 7-2.0  ft  during 
1978  to  1984. 

53.  Based  on  results  from  the  LAEMSED  numerical  model,  there  would  be 
about  a  20  percent  decrease  in  the  shoaling  rate  associated  with  removal  of 
the  six  major  withdrawals  and  discharges  in  the  harbor. 

54.  Advance  maintenance  of  6  ft  during  dredging  in  the  harbor  entrance 
and  Corpus  Christ!  Turning  Basin  would  increase  shoaling  rates  by  a  maximum  of 
about  1  in.  per  year.  In  the  areas  with  the  greatest  shoaling  rate,  no  in¬ 
crease  in  shoaling  was  observed;  therefore,  advance  maintenance  would  lessen 
the  required  dredging  frequency.  An  addition  of  a  sill  in  combination  with 
advance  maintenance  did  not  decrease  sedimentation  in  the  inner  harbor  by  a 
significant  amount. 

55.  Results  from  the  Phase  I  reconnaissance  survey  and  the  Phase  II 
data  collection  appear  to  indicate  that  a  primary  source  of  shoaling  is  the 
dredged  material  disposal  sites  south  of  the  channel  closest  to  the  harbor 
entrance.  Field  data  and  the  TABS-2  hydrodynamic  model  have  shown  that  cur¬ 
rents  generally  flow  from  the  disposal  areas  toward  the  channel  most  of  the 
time.  Due  to  the  apparently  high  erodibility  of  the  dredged  material,  it  does 
not  remain  at  the  disposal  site,  but  flows  back  into  the  channel  and  Is  pulled 
into  the  harbor  by  net  inflow  currents  along  the  bottom.  A  disposal  monitor¬ 
ing  program  is  necessary  to  determine  the  rates  of  erosion  of  the  material  at 
these  sites.  It  is  recommended  that  periodic  surveys  of  the  disposal  sites  be 
performed,  especially  following  a  dredging  operation,  to  determine  the  rate  of 
erosion  and/or  consolidation  of  the  dredged  material. 

56.  Tests  using  the  TABS-2  models  have  shown  that  relocation  of  the 
disposal  sites  to  6,000  ft  north  of  the  channel  and  moving  them  slightly  east¬ 
ward  reduces  the  amount  of  sediment  entering  the  bay  channel  by  75  percent. 

All  indications  suggest  that  moving  the  disposal  areas  to  north  of  the  channel 
and  limiting  agitation  of  sediments  (e.g.,  by  shrimping  activity)  in  the  area 
will  reduce  the  amount  of  material  reentering  the  channel  and  thus  entering 
the  inner  harbor. 

57.  Model  results  and  field  data  analyses  (Appendix  A)  suggest  that  the 
shoaling  mechanism  consists  of  the  following  processes:  Bay  sediments. 


including  resuspended  dredged  material,  are  resuspended  and  travel  across  the 
channel  where  they  settle  into  the  deeper  water.  Circulations  in  the  channel, 
driven  by  tides,  density  currents  (due  to  hypersaline  water,  sediment  con¬ 
centration,  and  possibly  freshwater  discharges),  and  harbor  discharges  and 
withdrawals  draw  sediment-laden  water  into  the  harbor,  where  the  sediment 
deposits. 

58.  The  model  results  indirectly  suggest  that  one  of  the  following 
untested  remedies  might  reduce  sedimentation  in  the  harbor: 

a.  Protection  of  the  dredged  material  disposal  area  from  wave 
action  and  tidal  currents  by  construction  of  sills  along  two 
sides. 

b.  Trapping  of  sediment  by  a  dredged  trap  adjacent  to  the  harbor 
charnel . 

£.  Disruption  of  the  density  current  by 

(1)  A  physical  barrier  (flap  or  curtain). 

(2)  A  water  or  bubble  curtain. 

d.  Mechanical  dredging  instead  of  hydraulic  dredging  of  harbor. 

e.  Modification  of  channel  geometry  at  bridge  to  alter  deposition 
patterns. 

This  evaluation  does  not  consider  the  economics  of  such  solutions,  but  poses 
them  only  as  possible  technical  solutions. 

59.  Due  to  the  limited  time  available  to  conduct  the  study,  the  impact 
of  wind  fields  was  not  evaluated.  Additional  tests  to  evaluate  the  impact  of 
various  wind  conditions  on  the  alternative  plans  should  be  conducted  during 
the  detailed  design  of  improvements. 
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Coef  f icient 
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LAEMSED 


Critical  shear  stress  for  deposition, 
N/sq  m 

Critical  shear  stress  for  particle 
erosion,  N/sq  m 

Dry  weight  density-f reshly  deposited 
layer,  kg/cu  m 

Particle  specific  gravity 
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Table  2 

Industrial  Activity  Near  Harbor  Entrance 
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Amount 

mgd 

Discharge  or 
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Location  on 
Figure  2 
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(Municipal  Waste) 
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4 
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1 
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4 

Discharge 

2 

Southwestern  Refining 

1 
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3 
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Chemical 

I 
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4 

Central  Power  and  Eight 

450 

Withdrawal 

5 
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Chemical 

4 

Wit hdrawa 1 

6 

! 


0.0001 

0.0005 

* 

2.65 
0.00015 
1  5 


0.0001 

* 
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APPENDIX  A:  REPORT  ON  CONDITIONS  AFFECTING  SHOAI.ING ,  CORPUS 
CHRISTI  HARBOR,  TEXAS,  RECONNAISSANCE  SURVEY,  AUGUST  1984 


Introduction 


1.  Shoaling  rates  between  the  Corpus  Christi  Harbor  entrance  and  the 
Highway  181  bridge  and  in  the  turning  basin  of  the  harbor  have  been  reported 
to  be  360,000  and  550,000  cu  yd  annually,  respectively.  A  reconnaissance  sur¬ 
vey  was  conducted  in  August  1984  to  develop  information  on  the  causes  of  this 
shoaling,  the  characteristics  and  sources  of  the  sediments,  and  circulations 
in  the  harbor  and  Corpus  Christi  Bay. 

2.  A  previous  study  of  this  sedimentation  problem  was  performed  by  the 
Committee  on  Tidal  Hydraulics  (CTH) ,  Corps  of  Engineers,  US  Army  (1965).* 
Figure  A1  is  a  vicinity  sketch  of  the  harbor.  That  study  concentrated  on 
municipal/industrial  withdrawals  and  discharges,  and  on  chloride  conditions  in 
the  harbor.  The  findings  of  that  study  were  that  large  volumes  of  industrial 
withdrawals  and  relatively  fresh  discharges  created  a  circulation  between  the 
harbor  and  the  bay.  The  recommendations  were  to  further  define  water  usage  in 
the  harbor,  to  study  circulation,  and  to  measure  suspended  sediment  loads. 


Figure  A1 . 


Corpus  Christi  harbor 


sus- 


Ob  ject Ives 

3.  This  report  describes  the  results  of  salinity,  temperature, 
pended  sediment  concentration,  settling  velocitv,  and  bed  material  property 
measurements.  They  are  intended  to  more  clearly  define  the  sources,  trans¬ 
port,  and  deposition  of  sediments  in  the  harbor. 

Approach 


A.  Salinity,  temperature,  and  suspended  sediment  were  measured  to 
understand  circulations  and  mixing.  Suspended  sediments  and  currents  were 
used  to  estimate  transport  and  settling  rates.  Since  settling  is  such  an 
important  property,  it  was  also  determined  from  both  field  and  laboratory 
tests.  Bed  material  samples  were  taken  from  the  bay  to  identify  the  nature  of 
shoaling  material  and  possible  sources  of  material.  Tests  on  these  samples 
provided  information  on  sediment  characteristics  during  sedimentation  and 
hindered-sett ling  consolidation. 

Description  of  Tests 
Field  procedures  and  equipment 

5.  The  survey  took  place  20-23  August  1984.  Conditions  for  the  survey 
were  about  normal  for  the  time  of  year.  The  weather  was  warm  and  drv.  After¬ 
noon  winds  were  typically  southeast  12-15  knots.  Early  morning  winds  were 
south  3-5  knots.  Because  of  a  prolonged  dry  spell,  river  flows  were  verv  low. 

6.  A  21-ft  survey  boat  was  used.  Strong  winds  from  the  southeast  made 
working  in  the  bay  in  the  afternoon  impractical.  Sampling,  therefore,  was 
done  in  the  bay  in  the  early  morning  and  shifted  to  the  harbor  by  afternoon. 

7.  A  tide  gage  was  established  at  the  US  Army  Engineer  District, 
Galveston,  dock  (Engineers  dock,  Figure  All  near  the  entrance  to  the  harbor. 

A  Fischer  and  Porter  Model  1550  with  a  paper  punch  recorder  was  used.  It  was 
fitted  with  a  stilling  well  and  hydraulic  damper  to  exclude  short  waves.  No 
datum  tiane  was  established.  Water-surface  elevations  were  recorded  every 
15  min. 

8.  Two  current  meters  were  located  off  the  Engineers  lock.  One  was 
deployed  at  a  depth  of  6  ft  below  the  water  surface  and  the  other  f>  ft  •  * .  -r 


the  bottom  on  a  single  string.  The  water  depth  was  ?4  ft.  The  deplovment 
lasted  about  a  day.  F.ndeco  Model  174  recording  current  meters  with  conduct  tv 
itv  and  temperature  sensors  were  used.  Currents  were  recorded  at  ?-min  inter 
vals.  Threshold  speed  for  these  meters  is  0.08  ft /sec  and  the  manufacturer's 
stated  accuracy  is  13  percent  of  full  scale  (8  ft /seel.  Directional  accuracy 
is  17.1  deg  when  velocities  are  above  threshold. 

8.  Water  column  measurements  were  repeatedly  taken  at  a  network  of  har 
hor  and  bay  sampling  stations.  Tahle  A1  gives  the  channel  mile  locations  and 
description  of  the  harbor /channe 1  sampling  stations  from  the  Viola  Turning 
Basin  through  the  harbor  and  into  Corpus  Christ i  Bay  channel  to  dav  marker  A? 
Table  A2  gives  a  description  of  the  hav  stations.  Figure  A 1’  shows  their 
approximate  locations  in  the  bav. 

10.  Conductivity,  salinity,  temperature,  turbidity,  and  depth  were  mea¬ 
sured  in  situ  using  an  InterOceans  Model  SI  3  CSTP-Tr  measuring  device.  The 
instrument  consists  of  a  sensor  package  and  an  onboard  readout  unit.  The 


manufacturer  gives  its  precision  as  follows:  conductivity  t0.02  mmhos/cm; 
salinity  i0.05  ppt;  temperature  i0.02°  C;  and  depth  +  0.15  m. 

11.  Suspended  sediment  samples  were  taken  with  either  a  pump  or  water 
bottle  samplers.  The  former  consists  of  a  1 2-v  d-c  pump  and  100  ft  of 
1/4-in. -ID  plastic  tubing.  The  horizontal  intake  points  into  the  flow. 

Pumped  samples  were  stored  in  8-oz  bottles.  A  total  of  222  samples  were 
collected. 

12.  The  water  bottle  sampler  is  a  long  polyvinyl  chloride  tube  with 
closures  at  both  ends.  It  is  <.'  nloyed  from  the  survey  boat  suspended  from  a 
wire  and  sinker  in  the  horizontal  position.  The  sampler  orients  itself  with 
the  flow.  A  messenger  weight  is  sent  down  the  wire  to  close  the  device.  This 
sampler  was  used  to  collect  large-volume  (5-gal)  samples  and  for  settling 
velocity  determinations  in  the  field. 

13.  Six  field  settling  tests  were  performed  using  the  water  bottle  sam¬ 
plers.  After  collection  of  a  sample,  the  sampler  was  placed  vertically  in  a 
rack.  Subsamples  were  taken  from  the  sampling  tube  at  0,  R,  15,  30,  50,  120, 
180,  and  240  min.  These  subsamples  were  analyzed  for  concentration  of  sus¬ 
pended  material  (CSM)  and  the  data  treated  as  laboratorv  settling  test 
results. 

14.  Current  velocities  were  also  measured  from  the  survey  boat  using 
over-the-side  equipment.  The  current  meter  was  a  Curlv  Model  655  which  uses  a 
vertical-axis  cup-tvpe  impeller.  Accuracy  is  *0.1  ft/sec  for  current  speeds 
less  than  1  ft/sec  and  4 5  percent  for  current  speeds  greater  than  1  ft 'sec. 
Accuracy  of  the  direction  indicator  is  M0  deg  at  current  speeds  greater  than 
0.5  ft /sec.  Accuracv  of  the  current  measurements  also  depends  on  wave  action 
and  the  motion  of  the  survev  vessel. 

15.  Current  measurements  were  made  while  the  survev  boat  was  tied  t o  a 
pile  or  other  object.  Attempts  were  made  to  use  anchors  to  hold  the  boat  sta¬ 
tionary,  but  these  were  unstjccessf u  1  . 

15.  Bed  sediment  samples  were  taken  wrfth  a  I’et  erson-t  vpe  crab  sampler. 
Attempts  were  also  made  to  use  a  gravity  corer.  These  a: tempts  'ailed  as  will 
he  discussed  later.  Twentv-six  bed  sediment  samples  were  collected.  Their 
1 o cat  ions  are  given  in  Tahle  Ad. 
laboratory  and  data  analyses 

17.  Water  samples  were  analyzed  for  total  suspended  mater  i  a and  sa  1  i  n  - 
itv  or  were  vised  in  settling  tests.  Determination  of  CSV  wa-  fv  fill  rat  ion 


and  gravimetric  analyses.  Nuclepore  polycarbonate  filters  with  0.45-u  pore 
size  were  used.  Samples  were  drawn  through  the  filters  with  a  vacuum  system. 
The  filters  and  holders  were  then  washed  with  distilled  water  to  remove  salts. 
The  filters  were  dried  at  105°  C'  for  1  hr.  CSM  was  calculated  from  the  weight 
of  material  retained  on  the  filters  and  the  sample  volume  filtered. 

18.  Water  column  data  consisting  of  CSM  and  current  velocity  with  depth 
were  fit  t  c>  exponential  and  logarithmic  curves,  respectively.  Depth-averaged 
CSM  and  current  speed  l’  ,  friction  velocity  l'*  ,  suspended  sediment  trans¬ 
port  rate  Q(sed)  ,  and  CSM  stratification  were  then  computed.  CSM  stratifi¬ 
cation  is  equal  to  the  surface-to-bottom  difference  divided  bv  the  depth- 
averaged  CSM  and  is  related  to  the  settling  characteristics  of  the  sediment 
and  the  turbulence  characteristics  of  the  flow. 

19.  Salinitv  and  temperature  data  were  used  to  compute  water  density 
using  Bennett's  function  (Bennett  1976).  The  density  parameter  c  was 
used  for  the  density  due  to  salinitv  and  temperature,  where  r  equals 
(density  -  1 )  *  1,000.  (If  o,.-  equals  24.5,  then  density  equals 

1.0745  g/ml,  for  instance.)  The  density  anomalv  due  to  the  presence  of  sus¬ 
pended  sediment  was  also  computed  assuming  the  density  of  the  sediment  parti¬ 
cles  to  be  2.65  g/ml.  It  was  designated  ”-C  and  was  used  to  evaluate  the 

sm 

effect  of  suspended  material  on  circulation.  The  total  density  parameter  was 

then  a  -C 

ST  sm 

20.  Settling  tests  were  performed  on  three  5-gal  samples  in  the  labora¬ 
tory.  These  samples  were  vigorously  shaken  for  5  min  bv  hand  prior  to  being 
introduced  into  a  4-in. -diatr  bv  b-ft-high  settling  column.  Samples  were  with¬ 
drawn  from  near  the  bottom  of  the  tube  with  time,  and  these  subsamples 

ana  1 vzed  f  or  CSM . 

21.  Settling  tests  were  also  performed  on  lied  material  resuspended  into 
the  supernatant  of  the  5-gal  samples .  The  pr< < eduro  used  was  the  same  as  for 
the  other  settling  tests.  Besuspended  samples  were  vigorously  mixed  and 
poured  into  the  settling  column.  Samples  were  withdrawn  with  t i me  and  ana¬ 
lyzed  (or  (SM.  The  sediment  used  tor  these  tests  was  a  mixture  from  three 
stations  <  sediment  I. notions  14,  15,  and  1  b  described  it'  Tab  1  o  A1'  bunted  in 


an  area  ranging  from  under  the  bridge  to  between  the  breakwite 
nine  tests  w.i>  conducted  with  t  he  initial  CSM  v.;r  ied  hetwi  en 
b  ,  ‘*00  mg'..  These  tests  de  t  e  r  mi  :  •  !  Mr  dcfe*'.  bun  e  -  e  t  r  !  ;  •  v 
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22.  Observations  were  made  in  the  laboratory  of  the  thickness  of 
deposited  material  with  time  during  settling  tests  of  initial  CSM  values  of 
1  , 2  b  0  mg/'  and  above.  These  observations  were  used  in  coniunct ion  with  the 
results  of  the  settling  tests  to  compute  the  initial  densitv  of  the  bed  formed 
by  sediment  settling  from  suspension. 

23.  Raw  data  from  field  and  laboratory  settling  tests  consisted  of  CSM 

at  sampling  times.  A  regression  equation  was  fit  to  the  data  by  the  least 
squares  method.  Coefficients  from  this  equation  were  then  used  to  calculate 
the  settling  velocity  distribution.  Tests  with  initial  CSM  of  2,230  mg f  >’  and 
above  exhibited  interference  among  particles  which  made  settling  velocities 
almost  uniform.  For  these  tests,  Owen's  method  of  visually  following 

the  descent  of  the  strong  concentration  gradient  was  used  to  estimate  settling 
v  e  1  o  c  i  t  v  . 


24.  Bed  sediment  samples  were  analyzed  for  bulk  wet  densitv  (RUT'', 
cation  exchange  capacity  < CF.O ,  loss  on  ignition,  percent  moisture,  and  grain 
size  as  well  as  settling  velocity  as  noted  earlier.  BVD  was  determined  using 
23-ml  wide -mouth  pycnometers.  CF.C  was  determined  (after  air-drvinc  samples'! 
bv  a  sodium,  saturation  method.  1  ess  on  ignition  was  calculated  from,  weight 
loss  after  beating  samples  for  1  hr  at  SSO°  C.  Percent  moisture  was  deter¬ 
mined  bv  evaporation  at  l OS"  C.  Salts  were  not  removed .  (.rain  size  was 
determined  on  dispersed  samples  using  hydrometer  and  sieve  analyses. 

25.  CFC  was  used  as  a  measure  of  sediment  robes! veness :  the  higher 
the  CFC,  the  more  cohesive  the  sediment.  Typical  values  in  meq  ■  1  OP  g  for 
minera  logi  ca  i  groups  are  as  follows:  kaolinite,  '  - 1  R ;  illite,  I0-41''; 
chlorite,  10-40;  and  mont mor i 1 1 oni t e ,  SO- ISO. 

/S.  Consolidation  or  florculant  settling  tests  were  performed  on  the 
same  sediment  used  in  the  settling  tests.  Sediments  were  mixed  with  ".at  ive 
water  fo  com  ent  rat  ions  ranging  from  20  to  200  g/i.  Thev  were  poured  into 
t  ransparent  cylinders  about  1.2  ft  high  with  diameters  ranging  : rom  2- i  'F  t  o 


u  in.  depending  on  the  test  concent  rat  ion.  At  these  high  concent,  rat  i  or.;-  .  t  he 
suspension  behaved  as  a  mass  and  col  lapsed .  A  sharp  interface  formed  bet  ween 
the  sediment  and  a  clear  overlving  laver.  The  rate  of  descent  of  this  inter- 

-  this  sediment  at  this 


fate  was  equal  to  the  hindered  settling  velocity  for 
eoncen t  ra t  i on . 
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Test  Results 


Water  column  measurements 

27.  A  plot  of  the  water  levels  at  the  Engineers  dock  is  shown  in  Fig¬ 
ure  A3.  Tides  were  diurnal  with  ranges  of  0.5  to  0.8  ft. 

28.  Vertical  profiles  of  concentration  C  and  current  speed  l’  are 

grouped  bv  location:  bay  and  interior  channel.  The  first  group  is  the  chan¬ 
nel  stations  with  results  presented  in  Plates  AI-A28.  The  velocity  measure¬ 
ments  for  these  profiles  are  suspect  or  were  not  attempted  at  all.  Tests 
results  for  stations  located  in  the  bay  are  shown  on  Plates  A29-A51 .  Profiles 

in  Plates  A1-A51  are  shown  in  terms  of  normalized  distance  Z/H  ,  where  Z  is 

the  height  above  the  bed  and  H  is  the  total  depth.  Most  of  these  profiles 
include  velocity  measurements.  A  method  being  developed  under  the  "Fine- 
Grained  Shoaling  in  Navigation  Channels"  work  unit  of  the  Improvement  of 
Operations  and  Maintenance  Techniques  research  program  was  applied  to  these 
data.  Stratification  of  CSM  and  shear  velocity  estimates  based  on  Manning's 
friction  assumption  were  used  to  calculate  settling  velocity.  The  mean  set¬ 
tling  velocity  for  25  data  sets  was  0.25  mm/sec. 

29.  Tables  A4-A7  contain  observed  salinity,  temperature,  CSM,  o  , 

-C  ,  and  -  -c  .  Each  table  contains  a  transect  either  along  the  harbor 
sm 
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channel  or  through  the  hay. 
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30.  Current  speed  is  shown  with  depth  for  the  bay  stations  in 

Plates  A29-A51.  Plots  of  the  direction  and  depth-averaged  speed  are  shown  in 

Figures  A4,  A5,  and  A6  for  21,  22,  and  23  August,  respectively.  Plots  showing 

depth-averaged  CSM  and  a„  -C  densities  for  22  August  are  presented  in  Fig- 

o  i  sin 

urea  A7  and  A8,  respectively. 

31.  Data  from  the  current  meter  string  at  the  Engineers  dock  indicated 
a  very  short  ebb  at  both  the  surface  and  bottom,  with  speeds  of  0.2  to 

0.3  ft/sec  lasting  only  about  1-2  hr.  Currents  during  the  flood  were  stronger 
(0.4  ft/sec)  at  the  bottom  than  at  the  surface  (0.3  ft/sec).  Duration  of  the 
flood  was  4-6  hr.  Since  the  condition  and  position  of  these  meters  were  un¬ 
certain  during  the  deployment,  detailed  data  are  not  presented. 

32.  Table  A8  contains  the  results  of  four  field  settling  tests.  Cumu¬ 
lative  settling  velocity  distribution  in  percent  and  a  differential  distribu¬ 
tion  by  settling  classes  are  shown.  Some  parameters  for  the  settling  velocity 
distribution  are  also  given.  The  weighted  average  of  the  distribution  is  the 
most  representative  of  these.  Geometric  mean  and  standard  deviation,  skew¬ 
ness,  and  kortosis  are  also  shown.  The  mean  average  settling  velocity  of 
these  four  tests  was  0.25  mm/sec. 

Laboratory  results 

33.  Table  A9  contains  the  results  of  the  laboratory  tests  on  three 
5-gal  field  samples.  These  samples  were  neither  kept  cool  nor  chemically 
preserved  prior  to  analysis,  and  the  results  should  be  viewed  with  some 
caution. 

34.  Results  of  the  laboratory  tests  on  resuspended  bed  material  are 
given  in  Table  A10  for  concentrations  ranging  from  250  to  1,840  mg/l.  On 
concentrations  between  2,250  and  8,900  mg/ l,  the  tests  were  conducted  by 
observing  the  rate  of  descent  of  the  strong  concentration  gradient  which 
occurs  as  settling  velocities  become  more  uniform.  The  results  of  those  tests 
are  presented  in  Table  All,  which  also  gives  the  initial  bed  density  results. 

35.  Results  of  settling  tests  on  the  dependence  of  settling  velocity  on 
CSM  showed  that,  compared  to  other  sediments,  shoaling  materials  exhibited  a 
rather  small  dependence.  This  could  be  because  these  materials  are  so  highly 
cohesive  (as  indicated  by  CEC)  that  they  form  tight,  rapidly  settling  aggre¬ 
gates,  even  at  low  concentrations.  The  range  of  settling  velocities  is  typi- 


The  standard  deviation  of  the  settling  velocity  distributions  is  gener- 
'  small,  indicating  that  aggregates  are  fairly  uniform. 


Figure  A6.  Depth-averaged  velocities,  23  August  1984 ,  flood  tide 


densities,  22  August  1984 


Figure  A8.  Depth-averaged  0ST-C 

36.  Initial  bed  densities  were  calculated  from  deposit  thicknesses  and 
sedimentation  rates.  The  deposits  grew  essentially  linearly  during  the  tests, 
reached  a  maximum  thickness,  and  then  thinned.  Table  All  shows  estimated  bed 
concentrations  and  densities  for  three  times:  when  50  percent  of  the  material 
had  settled  (when  available),  when  the  peak  bed  thickness  occurred,  and  at  the 
end  of  the  test. 

37.  Consolidation  test  results  are  summarized  in  Table  A12.  Initial 
and  final  concentrations  and  settling  rates  are  given.  Figure  A9  summarizes 
the  effect  of  concentration  on  settling  rates  from  settling  tube  and  consoli¬ 
dation  tests.  It  can  be  seen  that  the  rapid  falloff  in  hindered  settling 
velocity  occurs  between  0.1  and  0.2  g/cu  cm.  Hindered  settling  velocities  at 
newly  deposited  concentrations  are  about  0.025  mm/sec.  A  characteristic  time 
T  for  the  newly  deposited  material  to  reach  fully  settled  conditions  can  be 

approximated  as  T  =  8.5H  where  T  is  in  hours  and  H  is  the  initial 

o  o 

height  of  the  deposit  in  metres. 
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CONCENTRATION,  G/CU  CM 


Effect  of  concentration  on  settling  velocity  from  settling 
and  consolidation  tests  * 


Discussion  of  Results 

Critique  of  tests 

38.  The  basic  assumption  of  the  reconnaissance  survey  was  that  condi¬ 
tions  sampled  were  representative  of  conditions  responsible  for  the  shoaling. 
This  implies  that  the  shoaling  is  nearly  continuous  or  constant  and  not 
dominated  by  extreme  events  such  as  storms  or  freshets.  This  assumption  will 
be  tested  later  by  comparing  calculated  to  observed  shoaling  rates. 

39.  Current  speeds  measured  from  the  survey  boat  were  found  to  be  very 


low  in  the  harbor.  As  mentioned  previously,  attempts  to  make  measurements 
from  an  anchored  survey  boat  were  not  successful.  The  motion  of  the  survey 
boat  was  apparently  sufficient  to  contaminate  the  true  current  at  the  harbor 
stations.  Even  when  the  survey  boat  was  moored  tightly  to  piers,  etc.,  the 
instrument  readings  were  erratic  and  sometimes  unreasonable.  Currents  were 
apparently  at  or  below  the  threshold  of  the  current  meter.  Results,  there¬ 
fore,  are  more  qualitative  than  quantitative.  Current  measurements  in  the  bay 
were  not  affected  by  this  problem. 

40.  The  current  meter  string  was  struck  and  dragged,  presumably  by  a 
fishing  boat.  The  string  had  been  moved  to  shallower  water  when  retrieved. 
Directions  in  the  data  were  erratic,  and  therefore  these  data  are  also 
qualitative . 

41.  Because  of  the  numbers  of  fishing  boats  in  the  bay  and  the  rela¬ 
tively  short  periods  spent  in  the  bay  during  the  survey,  only  one  drogue 
deployment  was  carried  out. 

42.  The  settling  tests  relied  on  accurate  CSM  determination.  Low  CSM's 
made  two  of  the  field  settling  tests  unreliable. 

43.  Core  sampling  of  bed  material  was  unsuccessful.  A  gravity  corer  of 
100  lb  was  used  on  a  handline.  The  device  that  retains  the  sample  in  the 
corer  did  not  work  properly. 

44.  Table  A13  presents  a  summary  of  the  CEC,  loss  on  ignition,  percent 
moisture,  BWD,  and  clay-silt-sand  distribution.  Plates  A52-A77  present  dis¬ 
persed  particle  size  distribution  plots  for  the  samples. 

Sources  of  sediments 

45.  CTH  (1965)  found  that  the  Nueces  River  carried  an  average 
300,000  cu  yd  of  sediment  annually.  The  material  was  much  like  the  shoaling 
material  in  grain  size.  The  reconnaissance  survey  also  found  some  bed  mate¬ 
rial  in  the  bay  similar  to  the  shoaling  material  in  grain  size  and  other 
characteristics,  but  not  in  BWD.  Some  signs  of  bank  erosion  on  the  northern 
edge  of  Corpus  Christi  Bay  were  noted  during  the  survey.  It  appears  that 
there  are  multiple  sources  that  supply  fine-grained  material  to  the  bay. 

46.  Oso  Bay  Sewage  Treatment  Plant,  which  is  similar  in  treatment  level 
and  magnitude  to  the  effluent  discharge  to  the  harbor,  discharges  about 
330,000  kg  per  year  of  waste  solids.  Because  this  is  two  orders  of  magnitude 
lower  than  the  shoaling  rate,  waste  solids  do  not  seem  to  be  contributing  to 
shoaling  directly. 


47.  The  disposal  of  dredged  material  in  the  bay  could  have  acted  as  a 
source  of  shoaling  material.  Hydraulically  dredged  fine-grained  material  has 
a  concentration  of  roughly  0.075  g/ml.  This  material  will  assume  a  finite 
angle  of  repose  only  after  it  densities  to  about  0.15  g/ml.  Therefore  the 
disposed  material  will  spread  over  a  large  area  before  consolidating.  The 
large  surface  area  encourages  the  loss  of  material  to  erosion.  Material  which 
is  consolidated  from  dense  suspensions  takes  longer  to  reach  a  fully  settled 
state  than  sediments  which  are  settled  from  low  concentration  suspension. 

Even  after  5  weeks  of  consolidation,  sediments  used  in  consolidation  testing 
had  reached  concentrations  of  only  0.24  g/ml  as  opposed  to  their  in-place  con¬ 
centrations  of  0.36  g/ml.  The  shear  strengths  of  these  materials  are  very 
strongly  dependent  on  concentration.  Some  studies  have  used  a  power  law  with 
an  exponent  of  2.5.  Therefore,  disposed  dredged  material  probably  remains 
more  susceptible  to  erosion  than  the  shoaled  material  from  which  it  comes  for 
a  long  period  of  time. 

48.  The  highest  CSM's  found  during  the  survey  were  in  Corpus  Christi 
Channel  in  the  bay.  The  highest  depth-averaged  CSM  was  at  DM#62  and  the 
highest  point  sample  CSM  was  at  the  breakwaters.  These  high  CSM's  were 
associated  with  the  channel  rather  than  the  shallow  areas  of  the  bay.  The 
channel  is  deep  enough  that  wind-wave  action  did  not  disturb  bed  sediments. 

The  characteristic  settling  time  for  the  channel  is  about  16  hr  based  on  set¬ 
tling  velocities. 

49.  CSM's  in  the  shallow  parts  of  the  bay  were  moderate.  Although  the 
afternoon  wind  conditions  probably  increased  CSM's  somewhat,  it  is  unlikely 
that  CSM  in  the  bay  would  have  compared  to  those  in  the  channel.  Currents  in 
the  shallow  portions  of  the  bay  were  strong  enough  to  have  kept  most  of  the 
fine-grained  material  suspended.  As  mentioned  in  paragraph  47,  fine-grained 
bed  sediments  of  the  shallow  bay  area  are  partially  consolidated,  with  BWD's 
greater  than  1.2.  No  repositories  of  newly  deposited  fine-grained  material  as 
fluid  mud  were  discovered. 

50.  Fine-grained  materials  entering  the  bay  from  the  Nueces  system  or 
by  erosion  of  bed  material  could  have  been  concentrated  in  the  deep  channel 
areas  by  several  processes.  Tidal  currents  in  the  channel  were  observed  by 
drogues  to  be  aligned  with  the  channel.  Currents  in  the  shallow  areas  adja¬ 
cent  to  the  channel  flowed  obliquely  to  the  channel.  The  surface  layer  of  the 
channel,  therefore,  may  be  supplied  with  an  advective  source  of  suspended 
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material.  Sediments  which  settled  out  of  the  surface  layer  could  have  been 
trapped  in  the  channel.  Also  hypersaline  waters  formed  by  evaporation  would 
tend  to  sink  into  the  depths  of  the  channel  and  may  have  carried  sediments 
with  them.  Salinities  associated  with  high  CSM  in  the  bay  channel  were  the 
highest  found  on  the  survey. 

51.  Figure  A10  shows  a  three-dimensional  plot  of  CSM  along  the  channel 
and  with  depth.  It  shows  graphically  the  relative  concentrations  observed  in 
the  channel.  The  data  collected  during  the  reconnaissance  survey  indicate 
that  sediments  in  suspension  in  the  bay  channel  periodically  intrude  into  the 
outer  harbor.  These  suspended  sediments,  therefore,  form  the  source  with  the 
most  direct  transport  route  to  the  harbor. 

52.  Figures  A11-A14  show  concentration  probability  plots  for  all  chan¬ 
nel  stations,  inner  harbors,  outer  harbor,  and  bay  channel,  respectively.  The 
ordinate  is  an  exponential  probability  function.  The  distributions  can  be 
seen  to  be  generally  lognormal,  although  large  surf ace-to-bottom  differences 
may  have  distorted  the  distributions  for  the  outer  harbor  and  bay  channel. 

Note  the  CSM's  were  generally  lowest  in  the  inner  harbor,  between  Avery  Point 
and  Viola  turning  basins  shown  In  Figure  A14,  and  that  there  is  a  change  in 
scale  for  this  figure. 

Transport  of  sediment 

53.  Sediments  responsible  for  the  shoaling  in  Corpus  Christi  Harbor  are 
transported  by  circulations  in  the  harbor  and  in  the  bay. 

54.  The  circulation  in  the  harbor  is  surprisingly  strong.  CTH  (1965) 
indicated  that  while  the  tidal  prism  was  120  mgd  during  its  survey,  the  flow 
into  the  harbor  was  3,000  mgd.  The  38-ft  depth  was  dominated  by  flows  enter¬ 
ing  the  harbor  during  the  entire  tidal  cycle,  and  the  30-ft  depth  was  domi¬ 
nated  by  flows  entering  the  harbor  88  percent  of  the  time.  The  3-  and  12-ft 
depths  were  dominated  by  flows  exiting  the  harbor. 

55.  The  density-driven  flows  out  of  the  harbor  are  equal  to  the 
freshwater  (waste)  discharges  times  their  dilution.  The  vertical  profiles 
of  salinity  from  the  reconnaissance  survey  were  used  to  estimate  the  dilu¬ 
tion  of  freshwater  inflows  to  the  harbor.  Dilution  can  be  estimated  by 


(Sg  -  Sa)/So  -  Sg)  where  is  the  salinity  of  the  effluent  discharged  into 
the  harbor  (assumed  to  be  3  ppt)  ,  is  the  inflow  salinity  at  the  bottom  of 
the  channel,  and  Sq  is  the  outflow  salinity  at  the  surface.  Based  on  eight 
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Figure  A14.  Weibull  probability  plot  of  Corpus  Christi 

Bay-Channel  samples 

profiles  from  between  the  breakwaters,  the  Engineers  dock,  and  the  bridge,  the 
average  dilution  was  estimated  as  A3. 

56.  A  flow  budget  was  constructed  assuming  that  the  tide  floods  at  the 
bottom  and  ebbs  at  the  surface,  that  the  freshwater  discharges  leave  the  har¬ 
bor  at  the  surface,  and  that  withdrawals  are  made  up  of  bottom  water.  The 
assumptions  about  the  discharges  and  withdrawals  follow  from  the  dilution 
argument  given  earlier.  Therefore,  if  only  about  2  percent  of  the  harbor 
water  was  composed  of  fresh  discharges,  only  about  2  mgd  of  the  fresh  dis¬ 
charges  would  have  been  withdrawn  from  the  harbor.  Using  the  industrial  with¬ 
drawals  and  discharges  compiled  in  the  CTH  (1965)  report,  the  following  flow 
budget  was  constructed: 

_ Flows,  mgd _ 


Source 

Surface 

Bot  tom 

Tidal  prism 

-284 

284 

Dilution  water 

-645 

645 

Fresh  discharges 

-15 

— 

Withdrawals 

— 

80 

Total 

-Q44 

1  ,009 

'  •-  'AV-W --V-' 
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The  net  f 1 ow  into  the  harbor  was  assumed  to  be  65  mgd  as  indicated  in  the  CTH 
report.  The  tidal  prism  is  based  on  the  average  tide  range. 

57.  The  largest  component  of  this  budget  can  be  seen  to  be  the  dilution 
water  flow.  Its  magnitude  is  sensitive  to  the  estimated  dilution  factor 
applied  to  the  freshwater  discharges.  This  flow  was  induced  as  a  direct 
result  of  the  freshwater  discharges  into  the  harbor. 

58.  The  geometry  of  the  system  also  contributes  to  the  distribution  of 
net  flows  in  the  harbor.  The  channel  is  400  ft  wide  in  the  bay,  and  the 
breakwaters  are  630  ft  apart.  At  the  bridge  the  channel  narrows  to  300  ft  and 
then  opens  to  about  1,000  ft  in  the  turning  basin.  Constrictions  in  width 
have  been  found  to  produce  converging  near-bottom  net  currents  in  other  tidal 
flows,  and  have  been  associated  with  shoaling  areas.  The  shoaling  area  is 
centered  about  this  constriction,  with  little  shoaling  occurring  in  the  west¬ 
ern  end  of  the  turning  basin.  In  addition  to  producing  a  near-bottom  conver¬ 
gence  zone,  the  constriction  also  probably  enhances  mixing  in  the  harbor  and 
thus  increases  dilution  water  flows. 


59.  Circulation  in  the  bav  can  be  inferred  from  the  o„  -C  distribu- 
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tion  for  22  August  (Figure  A8) .  It  appears  that  circulation  is  radially  out¬ 
ward  from  the  harbor  entrance.  Isopleths  of  density  run  generally  north- 
south.  The  spot  measurements  of  currents  in  the  bay  are  for  flood  (21  and 
23  August)  (Figures  A4  and  A6 ,  respectively)  and  near  high  water  (22  August) 
(Figure  A5) .  They  indicate  that  typical  bay  current  speeds  are  on  the  order 
of  0.5  ft/sec  and  generally  counterclockwise  during  flood.  In  the  Rincon 
channel  area  (RCtfl-4,  Figure  A2) ,  currents  were  strongly  influenced  by  the 
Nueces  Bay  system.  Drogues  launched  at  CCBfll,  shown  in  Figure  A2,  adjacent  to 
the  channel,  closely  paralleled  the  channel  for  3.5  hr. 

60.  The  contribution  of  CSM  to  density-driven  flows  was  normally  verv 
small.  When  concentrations  reached  about  200  mg/ 5.,  however,  CSM  did  appear  to 
contribute  to  both  longitudinal  gradients  in  the  channel  and  to  vertical 
gradients.  Intrusion  of  bay-channel  waters  into  the  harbor  mav  be  enhanced  bv 
the  presence  in  the  bay  channel  of  high  concentrations  of  suspended  material. 
Deposition  of  sediment 

61.  Two  methods  were  used  to  estimate  shoaling  rates  in  the  harbor. 
Estimates  were  made  to  test  the  hypothesis  that  the  reconnaissance  survev 
represents  important  shoaling  condition: . 

62.  One  me  t  hod  (if  estimating  shoaling  usee  a  s,  d  i  . . .  '-udec  t  .  t  ’■(• 


approach  taken  in  the  CTH  (1965)  report.  The  CTH  approach  used  point  measure¬ 
ments  of  CSM  and  currents  at  a  single  station  on  the  eastern  edge  of  the  turn¬ 
ing  basin  over  a  cycle.  This  method  accounted  for  158,000  cu  yd  annual  shoal¬ 
ing  rate  as  opposed  to  a  rate  of  550,000  cu  yd  established  from  dredging 
records . 

63.  A  flow  budget  was  presented  previously.  When  that  budget  is  com¬ 
bined  with  eight  CSM  profiles  (three  flood,  three  ebb,  and  two  slack  water), 

the  mass  of  sediment  settling  can  be  estimated  by  (CSM*Q) .  -  (CSM*Q) 

in  out 

where  Q  is  the  inflow  is  at  the  bottom  and  the  outflow  at  the  surface.  With 
a  conversion  for  a  specific  weight  of  361  kg/cu  m  (22.5  lb  dry  weight/cu  ft) 
obtained  from  laboratory  analyses  of  BWD,  the  estimated  annual  shoaling  volume 
was  445,000  cu  yd  total  for  the  harbor.  The  total  harbor  shoaling  was  esti¬ 
mated  to  be  910,000  cu  vd  per  year  by  dredging  records  in  the  CTH  (1965) 
report . 

64.  The  second  estimate  of  shoaling  was  obtained  by  calculating  settling 
flux.  The  mass  rate  of  deposition  was  assumed  to  be  equal  to  the  bottom  CSM 
times  the  settling  velocity  (0.25  mm/sec).  It  was  assumed  that  all  sediment 
particles  which  settle  stick  to  the  bed. 

65.  For  the  outer  channel  between  the  breakwaters  and  the  bridge,  the 
average  near-bottom  concentration  was  99.5  mg/ J l.  The  estimated  average  shoal¬ 
ing  was  7.2  ft  per  year  or  an  average  shoaling  rate  of  1,000,000  cu  yd  per 
year. 

66.  The  turning  basin  had  an  average  near-bottom  CSM  of  28.2  mg /£.  This 
yielded  a  shoaling  thickness  of  2  ft  per  year  or  239,000  cu  yd  per  year.  The 
combined  shoaling  for  the  outer  channel  and  turning  basin  was  1,239,000  cu  yd 
per  year. 

67.  The  two  estimates  of  shoaling  rates  bracket  the  rates  reported  in 
the  CTH  report  (1965).  The  spread  of  these  estimates  is  not  unusual  and  is 
due  to  their  assumptions,  and  to  the  data  on  which  they  are  based. 

68.  The  conditions  found  during  the  reconnaissance  survey  appear  to  be 
representative  of  those  responsible  for  shoaling,  although  episodic  events  may 
also  be  important. 

Cone  1  us  ions 

69.  It  appears  that  the  shoaling  problem  is  caused  by  high 
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Table  A1 

Location  and  Nomenclature  of  Water  Column 
Sampling  Stations  In  the  Channel 


Channel 

Mile 

Location 

Abbreviation 

31.0 

Viola  Turning  Basin 

Viola 

29.2 

Tule  Lake  Turning  Basin  day  marker 

DM#  2 

26.7 

Tule  Lake  near  bulk  materials  dock  (day  marker) 

DM#  9 

24.9 

Avery  Point  Turning  Basin 

AP 

23.0 

Corpus  Christ!  Harbor  Turning  Basin 

TB 

22.6 

Highway  181  Bridge  overpass 

Brdg 

22.5 

Engineers  dock 

Eng  Dk 

22.2 

Between  breakwaters  at  harbor  entrance 

Jetty 

21.8 

Corpus  Christi  Channel  day  marker 

DM#82 

20.9 

Corpus  Christi  Channel  day  marker 

DM#  80 

19.6 

Corpus  Christi  Channel  day  marker 

DM#  7  4 

17.0 

Corpus  Christi  Channel  day  marker 

DM#  6  2 

Table  A2 

Location  and  Nomenclature  of  Water  Column 
Sampling  Stations  in  the  Bay 


_ Description _  Abbreviation 

Rincon  Channel 

Rincon  Channel  near  causeway  (pile)  RC/?1 

Rincon  Channel  day  marker  10  RC/?2 

Rincon  Channel  day  marker  6  RC/?3 

Halfway  between  RC/?3  and  DM#80  RC/?4 

Destroyed  pier  0.6  miles  north  of  breakwater  RC/?5 

Corpus  Christi  Bay 

Range  marker  south  of  DM/?80  CCB//1 

Sun  Oil  platform  CCB/?2 

Cities  Service  platform  CCB/?3 

Destroyed  range  marker  to  old  harbor  city  entrance  CCB/?4 

Gas  wellhead  CCB/?5 

Portland  CCB/?6 

Day  marker  south  and  east  of  Corpus  Christi  Small 

Boat  Harbor  CCB/?7 

Breakwater  opening  (day  marker  1)  CCB/?8 

Adjacent  to  Channel 


Day  marker 
Day  marker 
Day  marker 


DM/?  SOB 
DM#82B 
DM/?  86 


Table  A3 

Bed  Sample  Locations 


Sample 

No. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 


_ Location/Station _ 

Corpus  Christ!  Turning  Basin/TB 
Corpus  Christi  Turning  Basin 
Avery  Point  Turning  Basin/AP 
Chemical  Turning  Basin 
Tule  Lake  Turning  Basin 

Viola  Turning  Basin 

Viola  Turning  Basin/Viola 

Viola  Channel 

Tule  Lake  Channel 

Tule  Lake  Channel  -  North  Flank 

Tule  Lake  Channel  Bulk  Materials  Dock/DM#9 

Chemical  Turning  Basin 

Chemical  Turning  Basin  -  North  Flank 

Under  181  Bridge /Bridge 

Engineers  Dock 

Between  breakwaters/ jetties 

Rincon  Channel /RC01 

Rincon  Channel/RC02 

Rincon  Channel /RC# 3 

Between  Rincon  Channel  and  DM#80/RC#4 

Corpus  Christi  Channel  DM£80/DM#80B 

Corpus  Christi  Channel  DM#80  center  line/DM#80 

Corpus  Christi  Channel  DM#82/DM/(82B 

Corpus  Christi  Channel  center  line  at  DM#82/DM#82 

Between  breakwaters/ jetties 

Off  north  breakwater  at  DM#86 


Channel 

Mile 


22.8 

23.2 

24.9 

25.8 

29.1 

31.2 
31.0 

29.8 

27.3 

27.1 

26.7 

25.7 
25.7 
22.6 
22.5 

22.2 


20.9 

20.9 

21.8 

21.8 

22.2 


Table  A4 

Sample  Data  for  20  August  1984 


Time  Depth  Salinity  Temperature 


Station  CST 


ST  mg/ l 


o-C  ac„,-C 
sm  Sl  sm 


Time 

Depth 

Salinity 

Temperature 

CSM 

Station 

CST 

m 

ppt 

°C 

°ST 

mg/ £ 

a-C 

sm 

aST~Csm 

RC#1 

0615 

1.0 

37.8 

28.8 

24.3 

32 

0.0 

24.3 

0.3 

37.6 

28.9 

24.1 

33 

0.0 

24.1 

RC#2 

0635 

2.4 

37.6 

29.1 

24.0 

81 

0.1 

24.1 

0.6 

38.0 

28.9 

24.4 

65 

0.0 

24.4 

RC#3 

0645 

3.5 

37.4 

29.0 

23.9 

45 

0.0 

23.9 

2.1 

37.4 

29.0 

23.9 

45 

0.0 

23.9 

0.6 

37.6 

29.0 

24.0 

46 

0.0 

24.1 

RC#4 

0710 

3.5 

37.4 

29.0 

23.9 

41 

0.0 

23.9 

2.0 

37.4 

29.0 

23.9 

35 

0.0 

23.9 

0.6 

37.5 

28.9 

24.0 

19 

0.0 

24.0 

DM#80B 

0725 

4.0 

37.3 

29.1 

23.8 

41 

0.0 

23.8 

2.9 

37.3 

29.2 

23.8 

44 

0.0 

23.8 

0.6 

37.4 

29.1 

23.8 

28 

0.0 

23.9 

DM/?  80 

0730 

12.0 

37.3 

29.2 

23.8 

91 

0.1 

23.8 

7.1 

37.3 

29.3 

23.8 

55 

0.0 

23.8 

5.1 

37.3 

29.3 

23.7 

58 

0.0 

23.8 

0.6 

37.4 

29.2 

23.8 

33 

0.0 

23.8 

DM#82B 

0750 

4.1 

38.0 

29.1 

24.3 

96 

0.1 

24.3 

2.8 

37.3 

29.1 

23.8 

28 

0.0 

23.8 

0.6 

37.5 

29.3 

23.9 

31 

0.0 

23.9 

DM#  8  2 

0815 

13.7 

37.5 

29.1 

24.0 

62 

0.0 

24.0 

9.8 

37.7 

29.2 

24.0 

57 

0.0 

24.1 

5.6 

37.4 

29.2 

23.9 

42 

0.0 

23.9 

0.6 

37.6 

29.3 

23.9 

35 

0.0 

24.0 

Jetty 

0835 

13.5 

37.7 

29.2 

24.0 

44 

0.0 

24.1 

9.4 

37.6 

29.3 

24.0 

35 

0.0 

24.0 

5.0 

37.2 

29.4 

23.6 

54 

0.0 

23.7 

0.6 

36.9 

29.5 

23.4 

10 

0.0 

23.4 

Brdg 

0920 

13.0 

37.7 

29.3 

24.0 

41 

0.0 

24.0 

9.9 

37.6 

29.3 

24.0 

36 

0.0 

24.0 

5.4 

37.0 

29.6 

23.4 

16 

0.0 

23.4 

0.6 

36.8 

29.6 

23.3 

4 

0.0 

23.3 

(Cont inued) 
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Table  A5  (Concluded) 

Time  Depth  Salinity  Temperature 
Station  CST  m  ppt  °C 


Table  A6 

Sample  Data  for  22  August  1984 


Table  A7 

Sample  Data  for  23  August  1984 


A8  (Concluded) 
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Table  A12 

Consolidation  Rates 


Initial 

Settling  Rate 
mm/ sec 

Final 

Concentration, 
Dry,  g/cu  cm 

BWD 

g/cu  cm 

Time 

hr 

Concentration, 
Dry,  g/cu  cm 

BWD 

g/cu  cm 

0.0203 

1.0372 

0.1000 

2.1 

0.0738 

1.0700 

0.0290 

1.0423 

0.0508 

98.3 

0.1835 

1.1372 

0.0339 

1.0455 

0.0419 

98.1 

0.1803 

1.1353 

0.0460 

1.0529 

0.0328 

72.1 

0.1658 

1.1264 

0.0576 

1.060 

0.0284 

69.0 

0.1623 

1.1242 

0.0742 

1.0702 

0.0236 

68.4 

0.1746 

1.1318 

0.0804 

1.0740 

0.0229 

234.5 

0.2144 

1.1562 

0.1997 

1.1472 

0.00009 

168.0 

0.2336 

1.1680 
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APPENDIX  B:  THE  TABS-2  SYSTEM 


1.  TABS-2  is  a  collection  of  generalized  computer  programs  and  utility 
codes  integrated  into  a  numerical  modeling  system  for  studying  two-dimensional 
hydrodynamics,  sedimentation,  and  transport  problems  in  rivers,  reservoirs, 
bays,  and  estuaries.  A  schematic  representation  of  the  system  is  shown  in 
Figure  Bl.  It  can  be  used  either  as  a  stand-alone  solution  technique  or  as  a 
step  in  the  hybrid  modeling  approach.  The  basic  concept  is  to  calculate 
water-surface  elevations,  current  patterns,  sediment  erosion,  transport  and 
deposition,  the  resulting  bed  surface  elevations,  and  the  feedback  to  hydrau¬ 
lics.  Existing  and  proposed  geometry  can  be  analyzed  to  determine  the  impact 
on  sedimentation  of  project  designs  and  to  determine  the  impact  of  project  de¬ 
signs  on  salinity  and  on  the  stream  system.  The  system  is  described  in  detail 
by  Thomas  and  McAnally  (1985a).* 

2.  The  three  basic  components  of  the  system  are  as  follows: 

a.  "A  Two-Dimensional  Model  for  Free  Surface  Flows,"  RMA-2V. 

b.  "Sediment  Transport  in  Unsteady  2-Dimensional  Flows,  Horizontal 
Plane,"  STUDH. 

£.  "Two-Dimensional  Finite  Element  Program  for  Water  Quality," 
RMA-4. 

3.  RMA-2V  is  a  finite  element  solution  of  the  Reynolds  form  of  the 
Navier-Stokes  equations  for  turbulent  flows.  Friction  is  calculated  with 
Manning's  equation  and  eddy  viscosity  coefficients  are  used  to  define  the 
turbulent  losses.  A  velocity  form  of  the  basic  equation  is  used  with  side 
boundaries  treated  as  either  slip  or  static.  The  model  automatically  recog¬ 
nizes  dry  elements  and  corrects  the  mesh  accordingly.  Boundary  conditions  may 
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*  References  in  this  Appendix  are  located  in  the  References  section  of  the 
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be  water-surface  elevations,  velocities,  or  discharges  and  may  occur  inside 
the  mesh  as  well  as  along  the  edges. 

4.  The  sedimentation  model,  STUDH ,  solves  the  convection-diffusion 
equation  with  bed  source  terms.  These  terms  are  structured  for  either  sand  or 
cohesive  sediments.  The  Ackers-White  (1973)  procedure  is  used  to  calculate  a 
sediment  transport  potential  for  the  sands  from  which  the  actual  transport  is 
calculated  based  on  availability.  Clay  erosion  is  based  on  work  by  Ariathurai 
and  Partheniades  and  the  deposition  of  clay  utilizes  Krone's  equations 
(Ariathurai,  MacArthur,  and  Krone  1977).  Deposited  material  forms  layers,  as 
shown  in  Figure  B2,  and  bookkeeping  allows  up  to  10  layers  at  each  node  for 
maintaining  separate  material  types,  deposit  thickness,  and  age.  The  code 
uses  the  same  mesh  as  RMA-2V. 


5.  Salinity  calculations,  RMA-4 ,  are  made  with  a  form  of  the 
convective-diffusion  equation  which  has  general  source-sink  terms.  Up  to 
seven  conservative  substances  or  substances  requiring  a  decay  term  can  be 
routed.  The  code  uses  the  same  mesh  as  RMA-2V. 

6.  Each  of  these  generalized  computer  codes  can  be  used  as  a  stand¬ 
alone  program,  but  to  facilitate  the  preparation  of  input  data  and  to  aid  in 
analyzing  results,  a  family  of  utility  programs  was  developed  for  the  ;  oil  v 
ing  purposes: 

a.  Digitizing 

b.  Mesh  generation 

c.  Spatial  data  management 

d.  Graphical  output 

£.  Output  analysis 

f.  File  management 

g.  Interfaces 

h.  Job  control  language 
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dividing  the  area  of  interest  into  smaller  subareas,  which  are  called  ele¬ 
ments.  The  dependent  variables  (e.g.,  water-surface  elevations  and  sediment 
concentrations)  are  approximated  over  each  element  by  continuous  functions 
which  interpolate  in  terms  of  unknown  point  (node)  values  of  the  variables. 

An  error,  defined  as  the  deviation  of  the  approximation  solution  from  the  cor¬ 
rect  solution,  is  minimized.  Then,  when  boundary  conditions  are  imposed,  a 
set  of  solvable  simultaneous  equations  is  created.  The  solution  is  continuous 
over  the  area  of  interest. 

9.  In  one-dimensional  problems,  elements  are  line  segments.  In  two- 
dimensional  problems,  the  elements  are  polygons,  usually  either  triangles  or 
quadrilaterals.  Nodes  are  located  on  the  edges  of  elements  and  occasionally 
inside  the  elements.  The  interpolating  functions  may  be  linear  or  higher 
order  polynomials.  Figure  B2  illustrates  a  quadrilateral  element  with  eight 
nodes  and  a  linear  solution  surface. 

10.  Most  water  resource  applications  of  the  finite  element  method  use 
the  Galerkin  method  of  weighted  residuals  to  minimize  error.  In  this  method 
the  residual,  the  total  error  between  the  approximate  and  correct  solutions, 
is  weighted  by  a  function  that  is  identical  with  the  interpolating  function 
and  then  minimized.  Minimization  results  in  a  set  of  simultaneous  equations 
in  terms  of  nodal  values  of  the  dependent  variable  (e.g.  water-surface  eleva¬ 
tions  or  sediment  concentration) .  The  time  portion  of  time-dependent  problems 
can  be  solved  by  the  finite  element  method,  but  it  is  generally  more  efficient 
to  express  derivatives  with  respect  to  time  in  finite  difference  form. 

The  Hydrodynamic  Model,  RMA-2V 


11.  The  generalized  computer  program  RMA-2V  solves  the  depth-integrated 
equations  of  fluid  mass  and  momentum  conservation  in  two  horizontal  direc¬ 
tions.  .The  form  of  the  solved  equations  is 
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where 

h  =  depth 

u,v  =  velocities  in  the  Cartesian  directions 
x,t,y  =  Cartesian  coordinates  and  time 
p  =  density 

e  =  eddy  viscosity  coefficient,  for  xx  =  normal  direction  on 

x-axis  surface;  yy  *  normal  direction  on  y-axis  surface;  xy 
and  yx  =  shear  direction  on  each  surface 

g  =  acceleration  due  to  gravity 

a  «  elevation  of  bottom 

n  *  Manning's  n  value 

1.486  •  conversion  from  SI  (metric)  to  non-SI  units 
t  •  empirical  wind  shear  coefficient 
V&  -  wind  speed 
<|i  »  wind  direction 

w  *  rate  of  earth's  angular  rotation 
4>  *  local  latitude 

12.  Equations  B1 ,  B2,  and  B3  are  solved  by  the  finite  element  method 
using  Galerkin  weighted  residuals.  The  elements  may  be  either  quadrilaterals 
or  triangles  and  may  have  curved  (parabolic)  sides.  The  shape  functions  are 
quadratic  for  flow  and  linear  for  depth.  Integration  in  space  is  performed  by 
Gaussian  integration.  Derivatives  in  time  are  replaced  by  a  nonlinear  finite 
difference  approximation.  Variables  are  assumed  to  vary  over  each  time  inter¬ 
val  in  the  form 
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which  is  differentiated  with  respect  to  time,  and  cast  in  finite  difference 
form.  Letters  a  ,  b  ,  and  c  are  constants.  It  has  been  found  by  experi¬ 
ment  that  the  best  value  for  c  is  1.5  (Norton  and  King  1977). 

13.  The  solution  is  fully  implicit  and  the  set  of  simultaneous  equa¬ 
tions  is  solved  by  Newton-Raphson  iteration.  The  computer  code  executes  the 
solution  by  means  of  a  front-type  solver  that  assembles  a  portion  of  the 
matrix  and  solves  it  before  assembling  the  next  portion  of  the  matrix.  The 
front  solver's  efficiency  is  largely  independent  of  bandwidth  and  thus  does 
not  require  as  much  care  in  formation  of  the  computational  mesh  as  do  tradi¬ 
tional  solvers. 

14.  The  code  RMA-2V  is  based  on  the  earlier  version  RMA-2  (Norton  and 
King  1977)  but  differs  from  it  in  several  ways.  It  is  formulated  in  terms  of 
velocity  (v)  instead  of  unit  discharge  (vh) ,  which  improves  some  aspects  of 
the  code's  behavior;  it  permits  drying  and  wetting  of  areas  within  the  grid; 
and  it  permits  specification  of  turbulent  exchange  coefficients  in  directions 
other  than  along  the  x-  and  z-axes.  For  a  more  complete  description,  see 
Thomas  and  McAnally  (1985b). 

The  Sediment  Transport  Model,  STUDH 

15.  The  generalized  computer  program  STUDH  solves  the  depth- integrated 
convection-dispersion  equation  in  two  horizontal  dimensions  for  a  single  sedi¬ 
ment  constituent.  A  detailed  description  can  be  found  In  Thomas  and  McAnally 
(1985c).  The  form  of  the  solved  equation  is 
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where 


■  concentration  of  sediment 

“  depth-integrated  velocity  in  x-direction 

■  depth-integrated  velocity  in  y-direction 

*  dispersion  coefficient  in  x-direction 

*  dispersion  coefficient  in  y-direction 

«  coefficient  of  concentration-dependent  source/sink  term 

*  coefficient  of  source/sink  term 


16.  The  source/sink  terms  in  Equation  B5  are  computed  in  routines  that 
treat  the  interaction  of  the  flow  and  the  bed.  Separate  sections  of  the  code 
handle  computations  for  clay  bed  and  sand  bed  problems . 

Sand  transport 

17.  The  source /sink  terms  are  evaluated  by  first  computing  a  potential 
sand  transport  capacity  for  the  specified  flow  conditions,  comparing  that 
capacity  with  the  amount  of  sand  actually  being  transported,  and  then  eroding 
from  or  depositing  to  the  bed  at  a  rate  that  would  approach  the  equilibrium 
value  after  sufficient  elapsed  time. 

18.  The  potential  sand  transport  capacity  in  the  model  is  computed  by 
the  method  of  Ackers  and  White  (1973),  which  uses  a  transport  power  (work 
rate)  approach.  It  has  been  shown  to  provide  superior  results  for  transport 
under  steady-flow  conditions  (White,  Hilli,  and  Crabbe  1975)  and  for  combined 
waves  and  currents  (Swart  1976).  Flume  tests  at  the  US  Army  Engineer  Water¬ 
ways  Experiment  Station  have  shown  that  the  concept  is  valid  for  transport  by 
estuarine  currents. 

19.  The  total  load  transport  function  of  Ackers  and  White  is  based  upon 
a  dimensionless  grain  size 

v  -  °  (B6) 

where 

D  -  sediment  particle  diameter 
s  -  specific  gravity  of  the  sediment 
v  -  kinematic  viscosity  of  the  fluid 
and  a  sediment  mobility  parameter 
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where 
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total  boundary  shear  stress 

a  coefficient  expressing  the  relative  importance  of  bed-load  and 
suspended-load  transport,  given  in  Equation  B9 

boundary  surface  shear  stress 


The  surface  shear  stress  is  that  part  of  the  total  shear  stress  which  is  due 
to  the  rough  surface  of  the  bed  only,  i.e.,  not  including  that  part  due  to  bed 
forms  and  geometry.  It  therefore  corresponds  to  that  shear  stress  that  the 
flow  would  exert  on  a  plane  bed. 

20.  The  total  sediment  transport  is  expressed  as  an  effective 
concentration 
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where  U  is  the  average  flow  speed,  and  for  1  <  D  <_  60 
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For  D  <60 

gr 


n'  -  0.00 

A  -  0.17 

C  =  0.025 

m  =  1.5 
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21.  Equations  B6-B16  result  in  a  potential  sediment 
This  value  is  the  depth-averaged  concentration  of  sediment 
an  equilibrium  transport  rate  is  reached  with  a  nonlimited 
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where 


bed  shear  stress 


Cc  *  critical  concentration  ■  300  mg/Jt, 

If  the  bed  shear  stress  is  greater  than  the  critical  value  for  erosion  t 
material  is  removed  from  the  bed.  The  source  term  is  then  computed  by 
Ariathurai's  (Ariathurai,  MacArthur,  and  Krone  1977)  adaptation  of 
Partheniades'  (1962)  findings: 


h  (’. '  ‘) 


(B22) 


where  P  is  the  erosion  rate  constant,  unless  the  shear  stress  is  also 
greater  than  the  critical  value  for  mass  erosion  .  When  rg  is  exceeded, 
mass  failure  of  a  sediment  layer  occurs  and 
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where 

Tt  “  thickness  of  the  failed  layer 
P  *  density  of  the  failed  layer 

u 

23.  The  cohesive  sediment  bed  consists  of  1  to  10  layers,  each  with  a 
distinct  density  and  erosion  resistance.  The  layers  consolidate  with  over¬ 
burden  and  time. 

Bed  shear  stress 

24.  Bed  shear  stresses  are  calculated  from  the  flow  speed  according  to 
one  of  four  optional  equations:  the  smooth-wall  log  velocity  profile  or 
Manning  equation  for  flows  alone;  and  a  smooth  bed  or  rippled  bed  equation  for 


where 


f  =  shear  stress  coefficient  for  waves 
w 

u  =  maximum  orbital  velocity  of  waves 
om 

f  =  shear  stress  coefficient  for  currents 

A  Bijker-type  equation  for  total  shear  stress  caused  by  waves 
and  current 


1  ,  -1  7  1  f  2 

2  fcU  +  4  fwUom 


For  further  information  on  the  shear  stress  computation  equations,  see 
McAnally  and  Thomas.* 

Solution  method 

25.  Equation  B5  is  solved  by  the  finite  element  method  using  Galerkin 
weighted  residuals.  Like  RMA-2V,  which  uses  the  same  general  solution  tech¬ 
nique,  elements  are  quadrilateral  and  may  have  parabolic  sides.  Shape  func¬ 
tions  are  quadratic.  Integration  in  space  is  Gaussian.  Time-stepping  is 
performed  by  a  Crank-Nlcholson  approach  with  a  weighting  factor  (6)  of  0.66. 
A  front-type  solver  similar  to  that  in  RMA-2V  is  used. 


*  W.  H.  McAnally,  Jr.,  and  W.  A.  Thomas.  1980.  Shear  Stress  Computations 
in  Numerical  Model  for  Estuarine  Sediment  Transport,"  Memorandum  for  Record, 
US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 


APPENDIX  C:  THEORETICAL  ASPECTS  OF  LAEMSED 


1.  Under  the  Environmental  and  Water  Quality  Operational  Studies 
(EWQOS)  program  of  the  US  Army  Corps  of  Engineers,  developmental  work  was  con¬ 
ducted  on  a  two-dimensional  (2-D)  laterally  averaged,  f ree-surf ace ,  variable- 
density,  and  heat-conducting  model  for  use  in  simulating  flows  in  thermally 
stratified  reservoirs.  This  effort,  which  extended  the  earlier  work  funded  by 
the  US  Army  Engineer  Division,  Ohio  River  (Edinger  and  Buchak  1979) ,*  resulted 
in  a  numerically  efficient  model  that  is  known  as  LARM  (Laterally  Averaged 
Reservoir  Model).  Under  a  contract  with  the  US  Army  Engineer  District,  Savan¬ 
nah,  LARM  was  modified  for  use  in  computing  stratified  flows  in  estuaries  as  a 
result  of  both  salinity  and  thermal  effects.  This  model  is  known  as  LAEM 
(Laterally  Averaged  Estuarine  Model)  (Edinger  and  Buchak  1981) .  An  applica¬ 
tion  of  LAEM  to  a  channel  deepening  study  by  Johnson,  Boyd,  and  Keulegan 
(1987)  on  the  Lower  Mississippi  River  demonstrated  the  general  applicability 
of  LAEM  to  such  problems.  A  version  of  LAEM  which  allows  the  modeling  of 
interconnecting  channels  has  been  modified  by  Johnson  (Johnson,  Trawle,  and 
Kee,  in  preparation).  This  model,  which  computes  suspended  sediment  transport 
and  simulates  the  erosion/deposition  process  at  the  bed,  is  called  LAEMSED 
(Laterally  Averaged  Estuarine  Model  with  Sediment) . 

Governing  Flow-Transport  Equations 

2.  The  basic  set  of  flow  and  transport  equations  that  are  solved  in 
LAEMSED  are  statements  of  the  conservation  of  mass  and  momentum  of  the  flow 
field  plus  the  conservation  of  the  heat,  salt,  and  suspended  sediment  in  the 
water  body.  The  governing  equations  are  developed  by  first  performing  a 
temporal  averaging  of  the  three-dimensional  equations  for  laminar  flow. 
Boussinesq's  eddy  coefficient  concept  is  then  employed  to  account  for  the 
effect  of  turbulence  in  the  flow  field.  Next,  the  time-averaged  equations  are 
averaged  over  the  estuary  width  and  finally  over  an  individual  vertical  layer 
with  boundaries  at  z  =  k  +  1/2  and  z  =  k  -  1/2  where  z  is  the  positive 
downward  Cartesian  coordinate  and  k  is  the  integer  layer  number,  to  yield 


the  following  equations  that  are  solved  in  the  water  column  in  LAEMSED  (sym¬ 
bols  used  in  these  equations  are  defined  following  Equation  C9) : 


Heat  balance 
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Salinity  balance 
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Suspended  sediment  balance 
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Equation  of  state 
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-  5890  +  38T  -  0.375T  +  3S 
=  1779.5  +  11.25T  -  0.0745T2  -  (3.8  +  0.01T)S 


Variables  in  Equations  C1-C9  are  defined  as  follows: 
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x  and  z 


=  plan  area  of  overbank  storage,  m 

2 

=  x-direction  momentum  dispersion  coefficient,  m  /sec 

2 

=  z-direction  momentum  dispersion  coefficient,  m  /sec 
=  estuary  or  river  width,  m 

=  laterally  averaged  width  integrated  over  h  ,  m 

1/2 

=  Chezy  resistance  coefficient,  m  /sec 

=  resistance  coefficient  associated  with  wind 

3 

=  suspended  sediment  concentration,  kg/m 

2 

=  x-direction  temperature  and  salinity  dispersion  coefficient,  m  /sec 

2 

=  z-direction  temperature  and  salinity  dispersion  coefficient,  m  /sec 

2 

=  acceleration  due  to  gravity,  m/sec 
=  horizontal  layer  thickness,  m 

3 

=  source  strength  for  heat  balance,  °C  m  /sec 

3 

=  sediment  source,  kg/m  /sec 
=  integer  layer  number,  positive  downward 
=  pressure,  N/m^ 

3 

=  tributary  inflow  or  withdrawal,  m  /sec 

3 

=  exchange  of  flow  between  channel  and  overbank,  m  / sec 
=  laterally  averaged  salinity  integrated  over  h  ,  ppt 
=  laterally  averaged  salinity  in  top  layer,  ppt 
=  time,  sec 

°  laterally  averaged  temperature  integrated  over  h  ,  °C 
=  x-direction  laterally  averaged  velocity,  m/sec 

=  x-direction  laterally  averaged  velocity  integrated  over  h  ,  m/sec 

3 

=  cell  volume  (B  •  h  •  Ax) ,  m 
=  wind  speed,  m/sec 

=  z-direction  laterally  averaged  velocity,  m/sec 

=  Cartesian  coordinates:  x  is  along  the  estuary  center  line  at  the 
water  surface,  positive  to  the  right,  and  z  is  positive  downward 
from  the  x-axis ,  m 

=  longitudinal  spatial  step,  m 

=  time-step,  sec 

=  change  in  surface  elevation,  m 

=  specific  gravity  of  sediment 

=  surface  elevation,  m 
3 

=  density,  kg/m 


3 

p  =  air  density,  kg/m 
a 

t  -  water  density  times  tangential  stress  in  positive  x-direction, 

Z  2.  2 

m  f sec 

4>  =  wind  direction,  rad 

3.  Basic  assumptions  in  addition  to  the  reduced  dimensionality  are  that 
the  Boussinesq  approximation  (p  is  constant  except  where  multiplied  by  the 
acceleration  of  gravity)  is  applicable  and  that  vertical  accelerations  are 
negligible  so  that  the  pressure  can  be  considered  hydrostatic.  In  addition, 
the  concept  of  eddy  coefficients  is  used  to  represent  the  effect  of  both  time 
averaging,  as  previously  noted,  and  spatial  averaging  of  the  equations.  The 
horizontal  dispersion  coefficients,  A^  and  ,  are  assumed  to  be  constant, 
whereas  the  vertical  dispersion  coefficients,  and  ,  are  dependent  upon 

the  stratification  as  reflected  by  the  local  Richardson  number  ,  i.e., 


where 


£  i£ 


(Cll) 


and  A  and  D  are  the  vertical  coefficients  for  no  stratification.  Due 
z  z 

o  o 

to  the  hydrostatic  pressure  assumption,  unstable  stratification  cannot  be 

modeled  in  a  convective  fashion  and  thus  is  handled  in  a  diffusive  manner  by 

2 

increasing  D^_  to  its  stability  limit  of  h  /2At  ,  where  At  is  the  computa¬ 
tion  time-step. 

4.  The  laterally  averaged  horizontal  pressure  gradient  in  the  longitu¬ 
dinal  momentum  equation  contains  the  density  driving  force.  Using  the  expres¬ 
sion  for  the  hydrostatic  pressure,  the  horizontal  pressure  gradient  can  be 
divided  into  its  two  components  of  the  barotropic  (surface  slope)  gradient  and 
the  baroclinic  (density)  gradient  to  yield: 


to  -  -gp3x  "  g  *  ^  02 

Numerical  Solution  Scheme 


5.  Finite  difference  techniques  are  employed  to  solve  the  governing 
equations.  The  particular  scheme  employed  is  structured  such  that  the  water- 
surface  elevations  are  computed  implicitly.  Using  the  new  water-surface  ele¬ 
vations,  the  x-component  of  the  flow  velocity  is  then  explicitly  computed  from 
the  longitudinal  momentum  equation.  As  in  other  hydrostatic  models,  the  ver¬ 
tical  component  of  the  velocity  is  computed  from  the  continuity  equation  which 
is  reduced  to  the  incompressibility  condition  as  a  result  of  the  Boussinesq 
approximation.  The  solution  begins  at  the  bottom  and  progresses  up  the  column 
of  layers.  With  the  flow  field  computed,  the  water  temperature  and  salt  and 
suspended  sediment  concentrations  are  then  computed  from  their  respective 
transport  equations  in  a  semi-implicit  fashion.  Details  can  be  found  in 
Edinger  and  Buchak  (1979). 

6.  A  note  concerning  the  treatment  of  the  vertical  advection  term 

w  bC  |  -  Iw  bC  \  in  the  transport  equation  for  the  suspended  sedi- 

\  7  k+1/2  \  7  k-1/2 

ment  concentration  is  required.  In  the  equations  for  temperature  and  salin¬ 
ity,  the  vertical  advection  term  is  handled  explicitly.  However,  if  the 
settling  velocity  is  relatively  large  and/or  the  vertical  layer  spacing  is 
small,  an  explicit  representation  in  the  sediment  transport  equation  will  re¬ 
sult  in  a  severe  restriction  on  the  computational  time-step.  Therefore,  an 
approximate  factorization  scheme  has  been  used  to  implicitly  handle  the  verti¬ 
cal  advection  term  in  Equation  C8.  This  is  accomplished  in  the  following 
manner.  Before  averaging.  Equation  C8  had  the  form 

3C  3UC  3w,C 

3I1  +  IT  +  -3T2  -  R-H*s-  (cl3> 

where  R.H.S.  contains  all  other  terms  in  Equation  C8. 

Using 
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-  +  O(At) 


(Cl  4) 


at 


where  n  represents  the  time  level  and  O(At)  refers  to  additional  terms 
multiplied  by  At  and  higher  powers  of  At  and  substituting  into 
Equation  C13  yields 
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C°  1  +  At  +  — -  —  =  Cn  +  At  (R.H.S.)  +  0\At  / 

S  dX  o  2  S 


(C15) 


However,  the  left-hand  side,  written  in  operator  form,  can  be  factored  as 
(l  +  At  D’)  (l  +  At  D')  Cn+1  =  Cn+1  +  At  D*Cn+1  +  At  D'Cn+1  +  o(At2)  (C16) 

V  X/  \  2/  S  S  X  S  2  S 


where  and  are  differential  operators  defined  such  that 
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The  following  sequence  of  equations  can  then  be  written  for  Equation  C16. 


(1  +  At  Cg  =  c”  +  At  (R.H.S.) 


(C17a) 


(l  +  At  D’^  C"+1  -  Cg  +  o(At2) 


(Cl  7b) 


Thus,  with  Equation  C17a,  computations  are  first  made  along  the  channel  with¬ 
out  considering  the  vertical  advection  term.  This  solution  is  then  used  in  a 
vertical  sweep  with  Equation  C17b  to  yield  the  final  suspended  sediment  con¬ 
centration  field  at  the  new  time  level.  Unlike  standard  alternating- 
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direction-implicit  (ADI)  schemes,  to  an  accuracy  of  0(At  )  no  iteration  is 
required.  With  this  solution  procedure,  the  settling  velocity  of  the  sus¬ 
pended  sediment  does  not  influence  the  allowable  computational  time-step. 

7.  The  major  advantage  of  the  basic  solution  scheme  employed  in  LAEMSED 
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is  that  the  extremely  restrictive  stability  criterion  based  upon  the  speed  of 
the  free-surface  gravity  wave,  i.e, 


D 


where  H 


is  the 


At  <  (Ax/y-gTTax,  ......  max 

maximum  water  depth,  is  removed.  However,  since  the  convective  terms  and 
baroclinic  term  in  the  longitudinal  momentum  equation,  as  well  as  the  vertical 
diffusion  and  advection  terms  in  the  temperature  and  salt  transport  equations, 
are  lagged  in  an  explicit  fashion,  the  following  stability  criteria  still 
remain: 
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where  Ap  is  the  density  difference  of  the  fresh  and  saline  water,  and 
is  the  height  of  the  density  flow. 


w 


Governing  Bed  Equations 


8.  The  routines  in  LAEMSED  that  compute  the  exchange  of  sediment  be¬ 
tween  the  sediment  bed  and  the  water  column  are  modifications  of  those  found 
in  a  vertically  averaged  sediment  transport  model  called  STUDH  (Thomas  and 
McAnally  1985c).  The  sediment  may  be  treated  as  either  cohesive  (clay)  or 
noncohesive  (sand).  A  single,  effective  grain  size  is  considered  for  each. 
Conceptual  basis 

9.  The  following  are  basic  assumptions: 
a . 


Basic  processes  in  sedimentation  can  be  grouped  into  erosion, 
entrainment,  transportation,  and  deposition. 


b. 


Flowing  water  has  the  potential  to  erode,  entrain,  and  transport 
sediment  whether  or  not  sediment  particles  are  present. 
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c.  Sediment  on  the  streambed  will  remain  immobile  only  as  long  as 
the  energy  forces  in  the  flow  field  remain  less  than  the  criti¬ 
cal  shear  stress  threshold  for  erosion. 

d.  Even  when  sand  particles  become  mobile,  there  may  be  no  net 
change  in  the  surface  elevation  of  the  bed.  A  net  change  in 
the  surface  would  result  only  if  the  rate  of  erosion  was  dif¬ 
ferent  from  the  rate  of  deposition — two  processes  which  go  on 
continually  and  independently. 

e.  Cohesive  sediments  in  transport  will  remain  in  suspension  as 
long  as  the  bed  shear  stress  exceeds  the  critical  value  for 
deposition.  In  general,  simultaneous  deposition  and  erosion  of 
cohesive  sediments  do  not  occur. 

f.  The  structure  of  cohesive  sediment  beds  changes  with  time, 
overburden,  and  ambient  water  conditions. 

g.  The  major  portion  of  sediment  in  transport  can  be  characterized 
as  being  transported  in  suspension,  even  that  part  of  the  total 
load  that  is  transported  close  to  the  bed. 

Bed  shear  stress 

10.  Several  options  are  available  for  computing  the  shear  velocity  used 
in  the  bed  shear  stress  equation 

*  pu*  (C19) 

where 

p  «  water  density 
u^  -  shear  velocity 

a.  Smooth-wall  log  velocity  profile, 

/  u*h\ 

—  -  5.75  log  ^3.32  —J  (C20) 

where 

u  =  flow  velocity  in  bottom  computational  cell 
h  =  bottom  layer  thickness 
v  =  kinematic  viscosity  of  water 

b.  The  Manning  shear  stress  equation, 
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where 


n  -  Manning’s  roughness  value 


CME  «  coefficients  of  1  for  SI  (metric)  units  and  1.486 
for  non-SI  units 


c . 


A  Jonsson-type  equation  for  surface  shear  stress  (plane  beds) 
caused  by  waves  and  currents. 
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where 


f  *  shear  stress  coefficient  for  waves 
w 


u  «  maximum  orbital  velocity  of  waves 
om 


f  *  shear  stress  coefficient  for  currents 
c 


Bed  source  term 


11.  The  transport  equation  for  suspended  sediment,  Equation  C8, 


contains  a  bed  source  term,  H  .  The  form  of  this  term  is 

s 


a,  C  +  a., 

Is,  2 
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(C23) 


where  C 


Is  the  sediment  concentration  near  the  bottom  and  a.  and  a„ 
s,  1  2 

D 


are  computed  as  discussed  below.  Equation  C23  is  the  same  for  deposition  and 
erosion  of  both  sands  and  clays.  Naturally,  if  deposition  Is  occurring,  H 


will  be  negative,  whereas  a  positive  value  will  result  if  erosion  is  taking 


place.  In  Equation  C23,  is  in  units  of  1  per  second  while  a ^  is  the 


equilibrium  concentration  portion  of  the  source  term  and  is  in  units  of 


kilograms  per  cubic  metre  per  second.  Methods  of  computing  and 


depend  on  the  sediment  type  and  whether  erosion  or  deposition  is  occurring. 

12.  Sand  transport.  The  supply  of  sediment  from  the  bed  (i.e.,  the 
sediment  reservoir)  is  controlled  by  the  transport  potential  of  the  flow  and 
availability  of  material  in  the  bed.  The  bed  source  term  is 
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C  -  equilibrium  concentration 
t  »  characteristic  time  for  effecting  the  transition 
There  are  many  transport  relationships  for  calculating  C  for  saM.  The 
Ackers-White  (1973)  formula  was  adopted  for  STUDH  and  thus  is  also  used  in 
LAEMSED . 

13.  The  characteristic  time  t  is  somewhat  subjective.  It  should  be 

the  amount  of  time  required  for  the  concentration  in  the  bottom  layer  to 

change  from  C  to  C  .  In  the  case  of  deposition,  t  is  related  to 
s^  eq  c 

fall  velocity.  The  following  expression  is  employed: 

t  «  larger  of  C,  —  or  At  (C25 
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where 

*  coefficient  for  deposition 
W^  *  fall  velocity  of  sediment  particle 
At  ■  computation  time  interval 

In  the  case  of  erosion,  there  are  no  simple  parameters  to  employ.  The  follow 
ing  expression  is  used: 


t  -  larger  of  C  —  or  At 
c  e  u 


(C26 


where 

C  *  coefficient  for  entrainment 
e 

u  »  flow  velocity 

14.  Clay  transport.  Deposition  rates  of  clay  beds  are  calculated  with 


the  equations  of  Krone  (1962): 


where 


t  ■  bed  shear  stress 

=  critical  shear  stress  for  deposition 

C  *  critical  concentration  =  300  mg/ l 
c 

Erosion  rates  are  computed  by  Ariathurai's  simplification  (Ariathurai, 
MacArthur,  and  Krone  1977)  of  Parthenaides '  (1962)  results  for  particle  by 
particle  erosion.  The  source  term  is  computed  by 
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where 

P  =  erosion  rate  constant 

=  critical  shear  stress  for  particle  erosion 

15.  When  bed  shear  stress  is  high  enough  to  cause  mass  failure  of  a  bed 
layer,  the  erosion  source  term  is 
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where 

=  thickness  of  the  failed  layer 

Pj  =  density  of  the  failed  layer 
Atj.  *  time  interval  over  which  failure  occurs 
=  bulk  shear  strength  of  the  layer 

Bed  model 

16.  The  sink/source  term  in  Equation  C8  becomes  a  source/sink  term 

for  the  bed  model,  which  keeps  track  of  the  elevation,  composition,  and  char¬ 
acter  of  the  bed.  A  basic  assumption  is  that  the  width  over  which  deposition 
or  erosion  occurs  is  the  input  width  of  the  bottom  computational  cell. 

17.  Sand  beds.  Sand  beds  are  considered  to  consist  of  a  sediment 
reservoir  of  finite  thickness,  below  which  is  a  nonerodihle  surface.  Sediment 
is  added  to  or  removed  from  the  bed  at  a  rate  determined  by  a  weighted  value 
of  the  sink/source  term  at  the  previous  and  present  time-steps.  The  mass  rate 
of  exchange  with  the  bed  is  converted  to  a  volumetric  rate  of  change  by  the 
bed  porosity  parameter. 


Cl  2 


18.  Clay  beds.  Clay  beds  are  treated  as  a  succession  of  layers.  Each 
layer  has  its  own  characteristics  as  follows: 

a.  Thickness 

b.  Density 

c.  Age 

d.  Bulk  shear  strength 

e.  Type 

In  addition,  the  layer  type  specifies  a  second  list  of  characteristics: 

a.  Critical  shear  stress  for  erosion 

b.  Erosion  rate  constant 

£.  Initial  and  1-year  densities 

d.  Initial  and  1-year  bulk  shear  strengths 

e.  Consolidation  coefficient 

New  clay  deposits  form  layers  up  to  a  specified  initial  thickness  and  then 
increase  in  density  and  strength  with  increasing  overburden  pressure  and  age. 
Variation  with  overburden  occurs  by  increasing  the  layer  type  value  by  one  for 
each  additional  layer  deposited  above  it.  Change  with  time  is  governed  by  the 
equations 


jf(t0)  +  [f(tj) 
|f(tj)  +  M  log  l 


-  f(tn)]  log  (9t  +  1)  0  <  t  <  1  year 
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where 

f  -  time-varying  characteristics  of  density  or  bulk  strength 
t^  =  time  =  zero 
tj  *  time  -  1  year 
M  *  consolidation  coefficient 

Mass  deposition  rates  are  converted  to  volumetric  deposits  by  the  specified 
density  for  the  type  l  layer,  and  erosion  rates  are  converted  to  a  corre¬ 
sponding  volume  by  the  actual  density  of  the  eroding  layer. 
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Data  Requirements 


19.  The  major  data  input  required  by  LAEMSED  is  the  geometry  data 
describing  the  system.  At  the  center  of  each  computational  cell,  the  width  of 
the  estuary  or  river  must  be  prescribed.  In  addition,  overbank  areas  that 
contribute  to  the  storage  of  water  during  flood  tide  must  be  included  in  the 
model  geometry.  Other  data  required  are  the  boundary  conditions  that  drive 
the  internal  flow  field.  At  a  river  (upstream)  boundary,  a  discharge  hydro¬ 
graph  must  be  prescribed  along  with  the  temperature,  salinity,  and  concentra¬ 
tion  of  suspended  sediment  associated  with  the  inflow.  At  the  ocean  boundary, 
the  tide  must  be  prescribed.  At  tidal  boundaries,  vertical  distributions  of 
temperature,  salinity,  and  suspended  sediment  concentrations  must  also  be 
prescribed.  Water  temperature  is  often  assumed  constant  and  surface  heat 
exchange  is  usually  set  to  zero.  However,  for  problems  in  which  thermal 
effects  are  considered,  shortwave  solar  radiation,  air  temperature,  dew-point 
temperature,  and  wind  speed  must  be  known  in  order  to  compute  the  coefficients 
required  in  the  computation  of  the  rate  of  surface  heat  exchange. 

20.  In  addition  to  these  water  column  data  described,  information  about 
the  sediment  and  the  initial  bed  structure  must  be  input.  A  constant  settling 
velocity  is  currently  assumed.  Default  values  for  the  characteristic  param¬ 
eters  of  the  different  type  layers  that  can  make  up  a  sediment  bed  are  pro¬ 
vided  in  LAEMSED,  e.g.,  the  density  of  a  freshly  deposited  type  1  layer  is 

3 

defaulted  to  90  kg/m  .  However,  any  of  these  values  can  be  changed  through 
input  data,  if  desired. 
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